UNCLASSIFIED 


_ AD  NUMBER _ 

AD488395 

LIMITATION  CHANGES 
TO: 

Approved  for  public  release;  distribution  is 
unlimited. 


FROM: 

Distribution  authorized  to  U.S.  Gov't,  agencies 
and  their  contractors ; 

Administrative /Operational  Use;  JUN  1966.  Other 
requests  shall  be  referred  to  Office  of  Naval 
Research^  Arlington,  VA  22203. 


_ AUTHORITY 

ONR  Itr  27  Jul  1971 


THIS  PAGE  IS  UNCLASSIFIED 


COLUMBIA  UNIVERSITY 

22  o 

Ifeohnical  Report  No,  94 

REPRESENTATION  AND  REALIZATION  OF 
TIME -VARIABLE  LINEAR  SYSTEMS 

by 

Leonard  M.  Silverman 

ft?epared  for 

DEPARlMEWr  OF  THE  NAVY 
Office  of  Naval  Researohr 
Mathena^cs 

ItoMiington*  0.  c;  20'?60 

Cbn^ct  NpllR4^9(64) 

NR  044-ij.l 


June,  1966 


Reproduction  in  vAole  or  in  part  la  permitted  for  ar^ 
purpose  of  the  United  States  Government. 


DEPARTMENT  OE  ELECTRICAL  ENGINEERING 
SCHOOL  OF  ENGINEERING  AND  APPLIED  SCIENCE 
NEW  YORK  27,N.Y. 


f 


PRECIS 
Technical  Report 


Title:  "Representation  and  Realization  of  Time-Variable  Linear  Systems", 
Leonard  M.  Silverman,  Department  of  Electrical  Engineering  Technical 
Report  ^94,  Columbia  University;  NONR  4259(04). 

Background  :  This  report  is  an  outgrowth  of  o  study  of  nonlinear  ond  time- 
var table  systems  being  conducted  in  ihe  Department  of  Electrical  Engineer¬ 
ing,  Columbia  University.  Contained  in  the  report  are  the  results  of  an 
investigation  of  the  subclass  of  such  systems  described  by  a  system  of  first- 
order  linear  time-variable  differential  equations. 

Condensed  Report  Contents:  Systems  of  linear  time-variable  differential 
equations  are  studied,  with  particular  emphasis  on  identifying  those  system 
properties  and  concepts  that  can  be  characterized  without  knowledge  of 
the  equation's  solution.  Criteria  are  developed  for  determining  the  degree 
of  controllability  and  observability  of  such  systems,  These  criteria  are 
based  on  the  rank  of  matrices  formed  directly  from  the  system  coefficients. 
A  transformational  property  of  these  matrices  is  utilized  in  a  procedure 
for  reducing  noncontrol  lab  le  and  nonobservable  systems  to  lower  dynamic 
order.  Also  obtained  are  new  methods  for  characterizing  and  generalizing 
equivalent  system  representations,  including  criteria  for  equivalence  and 
zero-state  time  invariance  of  time-variable  systems. 

Based  on  the  theory  of  equivalent  systems,  a  new  approach  to 
the  synthesis  of  nonstationary  impulse  response  matrices  is  developed.  This 
method,  which  does  not  require  an  a  priori  assumption  of  separability,  pro¬ 
vides  a  systematic  procedure  for  reolizing  a  wide  class  of  responses. 

Application  is  also  made  to  time-variable  electrical  networks. 

By  relating  the  concepts  of  controllability  and  unilateral  transmission,  the 
existences  of  a  class  of  unilateral  networks  composed  solely  of  two-terminal 
RC  (time-variable)  components  is  established. 

For  Further  Information;  The  complete  report  is  available  from  the  Defense 
Documentation  Center  for  qualified  requestors. 
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ABSTRACT 


This  research  is  concerned  with  a  study  of  time -variable 
linear  systems  with  particular  emphasis  on  those  system  properties 
and  concepts  that  can  be  characterized  without  solving  time- variable 
differential  equations.  Criteria  are  developed  for  determining 
the  degree  of  controllability  and  observability  of  a  time- variable 
system.  These  criteria  are  based  on  the  rank  of  the  controllability 
and  observability  matrices,  newly  defined  quantities  formed  from 
the  system  coefficient  matrices  and  a  finite  number  of  their 
derivatives.  The  controllability  and  observability  matrices  are 
also  shown  to  be  useful  in  a  variety  of  other  system  analysis 
problems.  In  particular,  a  transformational  property  of  these 
matrices  is  utilized  in  an  explicit  method  for  reducing  noncontrollable 
and  nonobservable  systems  to  systems  of  lower  dynamic  order. 

Also  obtained  are  new  methods  for  characterizing  and  generating 
equivalent  system  representations,  including  criteria  for  equi¬ 
valence,  zero-state  equivalence  and  zero  state  time-invariance 
of  time-variable  systems.  Based  on  the  theory  of  equivalent  systems 
and  the  criteria  for  controllability  and  observability,  a  new  approach 
to  the  synthesis  of  nonstationary  impulse  response  matrices  is 
developed.  This  method  of  synthesis,  which  does  not  require  an  a 
priori  assumption  of  separability,  provides  a  systematic  procedure 
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for  realizing  a  wide  clasa  of  impulse  reaponses.  By  relating  the 
concepts  of  controllability  and  unilateral  transmission,  the  existence 
of  a  class  of  unilateral  networks  composed  solely  of  two-terminal 
RC  (time-variable)  components  is  demonstrated.  A  stable  example 
of  such  a  network  is  presented  and  possible  applications 
discussed. 
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CHAPTER  1 


INTRODUCTION 


1.  1  Motivation 

In  recent  years  there  has  been  a  growing  interest  in  the 
theory  and  application  of  systems  whose  parameters  vary  with  time. 
Among  the  types  of  systems  for  which  this  is  the  case  are  electrical 
networks,  control  systems  and  communication  channels.  Some 
examples  of  inherently  time >variable  networks  include  parametric 
amplifiers,  modulators  and  switched  networks.  In  addition,  problems 
involving  control  of  missile  flights,  or  communication  over  fading 
channels  often  may  be  analyzed  by  techniques  applicable  to  time- 
variable  systems.  A  fairly  general  survey  of  research  on  time- 
variable  systems  in  the  period  1950-1960  may  be  found  in  the  paper  of 
Zadeh  I  1  ].  For  contributions  prior  to  1950,  the  paper  of  Bennett [  2] 
is  valuable.  A  great  deal  of  emphasis  in  the  past  five  or  six  years 
has  been  on  those  systems  (both  fixed  and  time -variable)  that  can  be 
adequately  represented  in  state-variable  form[  3  ] .  This  has  especially 
been  true  in  modern  control  theory  (s^e  the  review  paper  of  Athans[4  ]), 
and  more  recently  in  network  analysis  [5-llJ. 

The  major  advantage  of  the  state-variable  approach  to  the 
analysis  of  time-variable  systems  is  that  an  essentially  geometric 
structure  is  given  to  the  differential  equations  which  describe  system 
behavior.  This  structure  considerably  reduces  notational  complexity 


and  furthers  insight  into  many  analysis  and  synthesis  problems. 

Although  there  is  a  fairly  complete  theory  of  time-variable 
systems  in  state- variable  form  [  3  ]  this  theory  is  for  the  most  part 
based  on  an  explicit  knowledge  of  system  solutions.  It  is  generally 
impossible,  however,  to  find  the  solution  to  a  set  of  time-variable 
differential  equations  in  closed  form.  While  it  is  possible  to 
compute  the  solution  to  any  desired  degree  of  accuracy  on  a  digital 
computer,  it  cannot  be  characterized  directly  in  terms  of  the 
system  parameters.  This  is  one  reason  why  relatively  little 
progress  has  been  made  in  the  analysis,  synthesis  and  application 
of  variable  parameter  systems  compared  with  what  has  been 
accomplished  for  fixed  systems. 

These  considerations  motivate  a  closer  examination  of 
several  important  properties  of  time-variable  systems  with  the 
intent  of  characterizing  them  as  completely  as  possible  in  terms  of 
’  own  system  parameters,  and  exploiting  them  in  several  specific 
applications . 

1 . 2  Problem  Formulation  and  Background 

The  major  system  concepts  examined  here  are  described 
below.  Precise  definitions  are  given  in  the  body  of  the  thesis. 

(1)  The  concept  of  controllability  is  essentially  concerned 
with  the  type  of  coupling  that  exists  between  input  and  state 
of  a  linear  system,  and  determines  the  extent  to  which  a 
system  can  be  controlled.  Dual  to  controllability,  the  concep'; 


of  observability  is  concerned  with  the  type  of  ..ouplinjy  *<^^1  < 
between  state  and  output  of  a  linear  systCTi,  and  deterrn; 
the  extent  to  which  a  system's  behavior  cr  a  be  obs'  -ved 

(2)  Related  io  controllability  and  oob'}tV^  'Uty 
the  idea  of  system  reducibility.  If  a  sy  scon  ’s  input- 
output  behavior  for  zero-initial  conditio  ib  cs-n  b'^  achieved 
with  a  system  of  lower  dynamic  order,  it  '3  said  i-,c  be 
reducible.  Noncontrollable  or  nonobservable  systems  are 
always  reducible. 

(3)  To  meaningfully  discuss  analysis  and  synthesis, 
some  notion  of  system  equivalence  is  essential.  Two 
systems  whose  states  are  related  by  a  nonsingular  trans¬ 
formation  of  coordinates  (possibly  time -variable)  are  said 
to  be  equivalent.  Also  important  is  the  concept  of  zero- 
state  equivalence.  Two  systems  having  the  same  input- 
output  response  for  zero  initial  conditions  are  said  to  be 
zero- state  equivalent.  Thus,  a  reduced  system  is  zero- 
state  equivalent  to  its  original  description. 

The  most  fundamental  of  the  above  system  properti)bs  are 
controllability  and  observability.  It  has  become  increasingly  evident 
in  recent  years  that  these  concepts  are  significant  in  diverse  problems 
of  system  analysis  and  synthesis  not  necessarily  related  to  control 
problems.  However,  controllability  first  played  an  important  role  in 


the  development  of  optimal  control  theory.  Criteria  which  may  be  viewed 
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as  types  of  controllability  were  present,  for  example  in  the  work  of 
Pontryagin  [iz],  Gamkrelidze  [  1  3  ] ,  LaSalle  [14],  and  Bertram  and 
Sarachik[  15],  Formal  definition  of  controllability  and  its  dual 
observability  was  made  by  Kalman  [  16-1 8  ],  who  also  showed  the 
importance  of  observability  in  optimal  filtering  problems  [l9,  20]  . 
Kreindler  and  Sarachik[2l  |  introduced  the  concept  of  output 
controllability,  and  indicated  some  subtleties  differentiating  various 
degrees  of  controllability  and  observability  in  the  time-variable  case. 
Unfortunately,  the  necessary  and  sufficient  conditions  formulated  by 
Kalman  and  his  colleagues  [  18,  19,  20,  22,  23]  ,  and  by  Kreindler  and 
Sarachik  [21  1  ,  for  the  various  forms  of  controllability  and  observability 
depend  explicitly  on  the  system  solution  in  the  time- variable  case. 

As  this  is  generally  not  available  except  as  a  numerically  tabulated 
solution,  it  has  been  difficult  to  analytically  characterize  the 
controllability  and  observability  properties  of  a  time-variable  system. 
The  same  is  true  for  reducibility  and  for  system  equivalence,  since 
known  conditions  for  these  properties  are  also  given  in  terms  of  system 
solutions. 

It  is  our  intent  to  obtain  criteria  for  controllability, 
observability  and  system  equivalence  which  in  no  way  depend  on 
knowledge  of  solutions  to  time-variable  differential  equations.  In 
addition,  methods  for  constructing  zero- state  equivalents  of  reducible 
systems  and  for  finding  explicit  transformations  between  equivalent 
systems  will  be  examined. 


Several  specific  problems,  whose  solutions  depend  on 
an  explicit  characterization  of  the  above  properties,  will  also  be 
attacked.  These  include: 

(1)  Conditions  for  the  existence  of  simple  canonical 
forms  for  time-variable  systems  and  methods  for  their  construction. 

(2)  A  systematic  procedure  for  realizing  a  prescribed 
impulse  response  matrix  as  a  system  of  differential  equations  in 
state-variable  form. 

(3)  Construction  of  unilateral  networks  from  two-terminal, 
time-variable  resistors,  inductors  and  capacitors. 

1 .  3  Summary 

In  Chapter  2  some  pertinent  information  from  the  theory  of 
linear  systems  is  reviewed  and  the  various  degrees  of  controllability 
and  observability  are  defined.  The  known  criteria  for  controllability 
and  observability  are  also  summarized  in  this  chapter. 

New  criteria  for  controllability  and  observability  are  derived 
in  Chapter  3.  These  criteria  are  based  on  the  rank  of  the 
"controllability"  and  "observability"  matrices,  which  are  formed  from 
the  system  coefficient  matrices  and  a  finite  number  of  their  derivatives. 
The  precise  extent  to  which  these  matrices  characterize  controllability 
and  observability  is  made  clear  by  relating  them  to  a  generalized 
Wronskian  matrix  of  vector  functions.  It  is  shown  that  this  matrix  has 
all  the  important  properties  possessed  by  the  common  Wronskian  of 
scalar  functions  [  24  ]  . 
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Reducibility  if  also  treated  in  Chapter  3.  The  controllability 
and  obiervability  matrices  are  used  not  only  to  test  for  reducibility 
but  in  the  actual  construction  of  reduced  equivalent  systems.  The 
methods  presented  for  system  reduction  are  shown  to  have  several 
important  advantages  in  comparison  with  previous  techniques  that  were 
restricted  to  either  fixed  systems  [  18]  or  a  special  class  of  time- 
variable  systems  [  25]  .  The  basis  of  the  present  reduction  scheme  is 
an  interesting  transformational  property  of  the  controllability  and 
observability  matrices.  This  property  is  quite  evident  once  observed, 
but  does  not  seem  to  have  been  exploited  previously,  even  for  fixed 
systems. 

F\irther  application  of  this  property  is  made  in  Chapter  4, 
where  a  general  approach  to  the  problems  of  system  equivalence  and 
representation  is  presented.  A  new  degree  of  controllability  and 
observability  is  defined  in  this  chapter.  For  the  class  of  systems 
possessing  these  properties  (which  includes  fixed  controllable  and 
observable  systems),  necessary  and  sufficient  conditions  for  system 
equivalence  are  derived,  together  with  explicit  methods  for  con¬ 
structing  transformations  relating  equivalent  systems.  Among  the 
interesting  and  potentially  useful  results  obtained  from  this 
investigation  is  a  criterion  for  determining  if  a  system  is  zero-state 
equivalent  to  a  fixed  system.  Also  given  in  Chapter  4  are  methods 
for  transforming  single-input- single-output  time-variable  systems 
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to  several  important  canonical  equivalents.  Included  are  the 
classical  input-output  differential  equation  form,  and  the  "phase- 
variable"  canonical  form  [  26-28  ],  widely  used  in  control  system 
applications.  These  canonical  structures  are  especially  valuable 
when  simulating  a  time- variable  system  on  an  analogue  (or  digital) 
computer  since  they  require  relatively  few  variable  components. 

The  theory  of  equivalent  systems  developed  in  Chapter  4 
provides  the  basis  for  a  new  approach  to  the  problem  of  synthesizing 
a  prescribed  impulse  response  matrix  as  a  system  of  differential 
equations  (i.  e.  ,  an  analogue  computer  realization).  This  method, 
presented  in  Chapter  5,  differs  significantly  from  previous  techniques 
[  18,  23,  29,  30]  in  that  an  a  priori  assumption  of  realizability  is 
not  required.  While  this  method  is  not  completely  general,  it 
provides  a  systematic  synthesis  procedure  for  a  wide  class  of 
responses.  In  particular,  for  stationary  impulse  response  matrices, 
a  theory  and  procedure  for  minimal  realization  is  obtained  which 
has  many  advantages  in  comparison  with  previous  approaches 
[  18,  30,  31,  32] . 

In  Chapter  6,  application  of  the  controllability  criteria  is 
made  to  a  problem  concerning  time-variable  networks.  While  it  is 
well  known  that  networks  containing  time-variable  RLC  components 
are  generally  non- reciprocal,  the  possibility  that  such  networks  may 
in  fact  exhibit  unilateral  behavior  does  not  appear  to  have  been  explored. 
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By  relating  the  lyetem  concept  of  controllability  to  the  network 
concept  of  unilateral  transmieeion,  we  are  able  to  eyatematically 
generate  a  clasi  of  unilateral  RC  networks. 
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CHAPTER  2 

TIME- VARIABLE  UNEAR  SYSTEMS 

2.  1  Basic  Definitions  and  Notation 

The  class  of  systems  to  be  considered  in  this  thesis  are 

those  describable  by  a  finite  set  of  first  order  differential  equations 
* 

of  the  form 


x(t) 

=  A(t)x(t)  +  B(t)u(t) 

(2.  la) 

y(t) 

=  C(t)x(t)  , 

(2.  lb) 

where  x(t),  an  n-vector,  is  the  state  of  the  system  at  time  t  ;  u(t), 
u.i  r-vector,  is  the  input:  and  y(t) ,  an  m-vector,  is  the  output.  The 
matrices  A(t),  B(t)  and  C(t),  are  of  order  compatible  with  the 
vectors  x(t),  u(t)  and  y(t). 


Notation:  Lower  case  symbols  will  be  used  to  denote  both  scalar 
and  vector  functions,  while  upper  case  symbols  will  be  reserved  for 
all  other  matrices.  When  the  context  is  clear,  the  explicit  dependence 
of  functions  on  their  argument  will  be  suppressed  (i.  e. ,  A  =  A(t)). 

The  operations  of  transposition  and  inversion  will  be  denoted  by  A' 
and  A~^,  respectively,  and  the  following  notations  will  be  used  for 
differentiation: 


•  *v 


r 


If  the  coefficient  matrices  are  sufficiently  well  behaved, 
the  output  of  such  a  system  is  given  by  [  1  ] 


i 

“•o'*  f 


y(t)  =  C(t)*(t,tQ)x(t^j)  +  I  C(t)*(t,r)B(r)u(r)dr  ,  (2.2) 


where  is  the  state  of  the  system  at  some  arbitrary  time  t^  , 

and  t(t,  r)  is  the  system  transition  matrix.  The  transition  matrix 
may  be  defined  as 


*(t,r)  =  X(t)X"  (T) 


(2.  3) 


where  X(t),  an  n  x  n  matrix,  is  a  fundamental  matrix  [  2  ]  of  (2.  la); 
that  is, 


X(t)  =  A(t)X(t) 


(2.4) 


and  X(t)  has  rank  n  for  all  t  . 


It  is  also  convenient  to  define  the  impulse  response  matrix 


of  system  (2.  1): 


r 

To  avoid  unnecessary  complication,  it  will  generally  be  assumed  that 
the  matrices  A,  B  and  C  are  continuous  functions  of  time. 
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H(t,  T)  = 


C(t)*(t,T)B(T), 


t  i  T 


t  <  0 


(2.  5) 


For  i.ro  initi.1  condttion.,  th.  impuloe  rerponre  matrix  completely 
characteriiee  the  input-output  behavior  of  iyetem  (2.  1)  eince  if 

xUq)  =  0  , 


y(t)  =  H(t,T)u(T)dT 


(2.6) 


From  (2.  3)  it  is  clear  that  the  impulse  response  matrix 


may  also  be  expressed  in  the  form 


H(t,T)  =  ’KOeiT)  , 


tax 


(2.7) 


where 


^(t)  =  C(t)X(t)  , 


(2.8) 


0(T)  =  X'  (T)B(T)  , 


(2.9) 


and  X(t)  i.  a  fundamenttl  matrix  of  (2.1a).  Although  d-ft)  and  e(t) 

are  not  unique  matricee  lor  a  given  eyetem,  they  are  unique  within 


a  constant 


tran.formation.  11  X,(t)  and  X^IO  are  any  two  lunda. 


r 


mental  matrices  of  (2.  la)  then 
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V., 

I 

h 


t 


Xj(t)  =  X^UIK 


(2.10) 


where  K  is  a  constant  nonsingular  matrix  [  2  ]  . 

It  is  often  usefiil  to  consider  various  equivalent  represent¬ 
ations  of  a  given  system  for  purposes  of  both  analysis  and  synthesis. 
There  are  many  types  of  system  equivalence  that  can  be  defi  id 
[  3,  4] ,  but  the  following  will  be  most  important  here. 

Definition  2.1;  Let  T(t)  be  an  n  x  n  matrix,  nonsingular  and 
continuously  differentiable  for  all  t  ,  and  let  z(t)  =  T(t)x(t)  . 

Then  it  will  be  said  that  the  system 


f  ■ 

ft. 


z  (t)  =  A(t)z(t)  +  B(t)u(t) 
y(t)  =  C(t)z(t) 


(2.11a) 

(2.  11b) 


.a-' 


where 


1 

A  =  (TA  +  f  )T‘^ 

(2.  12a) 

B  =  TB 

(2.  12b) 

t  . 

C  =  ct"^  , 

(2.  12c) 

i 


r-N-' 


1 

is  equivalent  (algebraically  equivalent  [3])  to  system  (2.  1)  and  that 
T  is  an  equivalence  transformation. 

If  X(t)  is  a  fundamental  matrix  of  (2.  la)  then  T(t)X(t) 
is  clearly  a  fundamental  matrix  of  (2.  11a),  so  that  the  transition 
matrix  of  (2.  11)  is  given  by 

I(t,T)  =  T(';$(t,T)T‘^T)  .  (2.13) 

The  impulse  response  matrix  is  thus  seen  to  be  invariant  under  an 
equivalence  transformation  since  from  (2.  12)  and  (2.  13) 

H(t,T)  =  C{t)J(t,T)BiT) 

=  C(t)$(t,T)B(r)  =  H(t,T)  . 

Furthermore,  if  x(tQ)  is  the  initial  state  of  (2.  1),  and  z(tQ)  =  T(tQ)x(tQ) 
is  the  initial  state  of  (2.  11)  then 

C(t)$(t,  tQ)z(tQ)  =  C(t)$(t,  tQ)x(tQ)  . 

Thus,  systems  equivalent  under  Definition  2.  1  are  both  zero-state 
and  zero- input  equivalent  [  4  ]  . 


'i7 
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2 ,  2  Controllability 

The  concept  of  controllability  as  it  arose  in  the  study  of 
optimal  control  problems  is  essentially  concerned  with  the  possible 
state  transitions  that  can  be  effected  in  a  system  by  application  of 
some  input.  There  are  two  basic  types  or  "degrees"  of  controllability 
having  importance  in  control  theory;  complete  [  3,  5,  6] ,  and  total 
[  7  ]  (differential  [  8  ]  ,  proper  [  9  ]  1  controllability.  The  major 
difference  between  the  two  is  that  total  controllability  insures  that 
the  state  of  the  system  can  be  controlled  as  quickly  as  desired, 
while  complete  controllability  only  implies  that  the  state  can  be 
controlled  in  some  finite  time. 

Precise  definitions  of  these  concepts  will  be  given  below. 
Various  equivalent  definitions,  as  well  as  a  fuller  discussion  of  the 
role  they  play  in  control  theory  may  be  found  in  the  work  of  Kalman 
[  3,  5  ]  ,  LaSalle  [  9  ]  ,  Kalman,  Ho  and  Narendra  [  6  ]  ,  Kreindler 
and  Sarachik  [  7  ]  ,  and  Weiss  and  Kalman  [  8  ]  . 

Definition  2.2  (Controllability): 

(a)  System  (2.  1)  is  said  to  be  completely  controllable  on  an  interval 
[  tQjtj]  if  for  any  state  at  t^  ,  and  any  desired  final  state 

at  tj  ,  there  exists  an  input  u(t)  defined  on  such  that 

x(tj)  =  Xj  . 

(b)  System  (2.  1)  is  said  to  be  completely  controllable  at  time  t^ 

if  there  exists  a  finite  time  t^  >  t  such  that  the  system  is  completely 


JTI-.  ^ 


'^rrtviww-.*  i 


controllable  on  [  tj]  •  system  (2.  1)  is  completely  controllable 
at  all  Iq,  it  is  said  to  be  completely  controllable. 

(c)  System  (2.  1)  is  said  to  be  totally  controllable  on  an  interval 
[  t^,  tj]  if  it  is  completely  controllable  on  every  subinterval  of 

[  tQ.tj]  . 

(d)  System  (2.  1)  is  said  to  be  totally  controllable  at  time  t^  if  for 
all  >  tp  ,  it  is  completely  controllable  on  [^Qjtj]  .  If  system  (2.  1) 
is  totally  controllable  at  all  t^  ,  it  is  said  to  be  totally  controllable. 

Necessary  and  sufficient  conditions  for  controllability  of 
a  time-variable  system  are  well  known  [  6,  7,  8]  and  will  be 
summarized  below.  If  X(t)  is  any  fundamental  matrix  for  (2.  la)  and 
<9(t)  =  X  ^(t)B(t),  then 

Theorem  2.1:  System  (2.1)  is  completely  controllable  on  the 

interval  [  t^,  t^]  if  and  only  if  the  rows  of  9(t)  are  linearly  inde¬ 
pendent  functions  of  T  on  [t^jt^]. 

Theorem  2.  2:  System  (2.  1)  is  totally  controllable  on  the  interval 

[  tp,  tj]  if  and  only  if  the  rows  of  9(T)  are  linearly  independent 
functions  of  T  on  every  subinterval  of  [  t^,  tj]  . 

Criteria  for  controllability  at  a  specific  initial  time,  and 
for  all  initial  times  follow  naturally  from  Definition  2.  1  and  the  above 
theorems  [7,8]. 

It  should  be  observed  that  it  is  usual  to  express  the 
controllability  criteria  in  terms  of  particular  fundamental  matrices 


■ '  *  *  ; 
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(either  'f(tj,T)  [7  ]  or  [  8  ]  ).  Since  any  two  fundamental 

matrices  are  related  by  a  constant  (with  respect  to  r  ),  nonsingular 
matrix  as  in  (2.  10),  the  independence  of  the  rows  of  0(T)  is  not 
affected  by  the  choice  of  fundamental  matrix  used  to  define  it. 

Furthermore,  controllability  is  invariant  under  an 
equivalence  transformation.  To  see  this,  let  X(t)  be  a  fundamental 
matrix  for  (2.  11a)  and  let 

0(t)  =  X‘Nt)B(t)  . 

By  definition,  B(t)  =  T(t)B(t),  and  as  noted  previously,  X(t)  =  T(t)X(t), 
where  X(t)  is  a  fundamental  matrix  for  (2.  la).  Therefore, 

0(t)  =  0(t)  and  anything  said  about  the  controllability  of  (2.  1),  must 
hold  for  (2.  11). 

2.  3  Observability: 

The  concept  of  observability  introduced  by  Kalman  [  3,  5,  10]  , 
is  of  considerable  importance  in  optimal  filtering  and  prediction  . 
problems.  Essentially,  if  a  system  is  observable,  its  state  at  some 
particular  time  can  be  determined  from  observations  of  the  system 
output.  Various  degrees  of  observability  will  be  defined  below.  It 
will  be  clear  from  these  definitions  and  the  criteria  which  follow  that 
observability  plays  a  role  dual  to  that  of  controllability  in  describing 
a  systems  structure.  This  duality  first  observed  by  Kalman  [  10  ] 


"T-- 
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greatly  simplifies  the  solution  of  many  problems  [  5,  10]  ,  Further 
discussion  of  observability  and  the  dual  relationship  it  bears  to 
controllability  may  be  found  in  [  7  ]  and  [  8  ]  . 

Definition  2.  3  (Observability): 

(a)  System  (2.  1)  is  said  to  be  completely  observable  on  an  interval 

[  tj]  initial  state  at  t^  can  be  determined  from 

knowledge  of  the  system's  output  over  [  t^,  t^]  . 

(b)  System  (2.  1)  is  said  to  be  completely  observable  at  time  t^  if 
there  exists  a  finite  time  tj  >  t^  such  that  the  system  is  completely 
observable  on  [  t^,  t^]  .  If  system  (2. 1)  is  completely  observable  at 
all  tg  ,  it  is  said  to  be  completely  observable. 

(c)  System  (2.  1)  is  said  to  be  totally  observable  on  an  interval 

[  tg,  tj]  if  it  is  completely  observable  on  every  subinterval  of  [  t^,  t^]  . 

(d)  System  (2.  1)  is  said  to  be  totally  observable  at  time  t^  if  for 
all  tj  >  t^,  it  is  completely  observable  on  [  t^,  t^]  .  If  system  (2.  1) 
is  totally  observable  at  all  t^  ,  it  is  said  to  be  totally  observable. 

If  X(t)  is  any  fundamental  matrix  for  (2.  la)  and 
^(t)  =  C(t)X(t)  then  the  following  conditions  for  observability  may  be 
established  [  7,  8  ]  . 

Theorem  2.  3:  System  (2.  1)  is  completely  observable  on  the  interval 
[  tQ,  tj]  if  and  only  if  the  columns  of  ’J^(t)  are  linearly  independent 
functions  of  t  on  [t^jt^]  . 


f 
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Theorem  2.  4;  System  (2.  1)  is  totally  observable  on  the  interval 
[  t^,  t^]  if  and  only  if  the  columns  of  ^'(t)  are  linearly  independent 
functions  of  t  on  ever  subinterval  of  [t^^tj]. 

The  form  of  Theorems  2.  3  and  2.  4  when  compared  with 
that  of  Theorems  2.  1  and  2.  2  immediately  suggests  a  precise 
formulation  of  the  duality  between  observability  and  controllability 

[7]. 

Theorem  2.  5:  System  (2.  1)  is  completely  controllable  (observable) 
on  the  interval  [  t^,  tj]  if  and  only  if  the  system 


w(t)  = 

-A'(t)w(t)  +  C'(t)u(t) 

(2.  14a) 

II 

B'(t)w(t) 

(2.  14b) 

is  completely  observable  (controllable)  on  [  t^,  t^] . 

To  prove  the  theorem,  it  is  only  necessary  to  observe 
that  if  X(t)  is  a  fundamental  matrix  of  (2.  la),  then  (X  \t))'  is  a 
fundamental  matrix  of  (2. 14a)  [  2  ] . 

If  system  (2. 1)  is  time- invariant,  all  degrees  of 
controllability  (and  observability)are  equivalent  and  we  may  classify 
systems  as  being  either  controllable  (observable)  or  non  controllable 
(non  observable) .  Most  importantly,  the  necessary  and  sufficient 
conditions  for  controllability  and  observability  in  the  fixed  case  can 
be  formulated  directly  in  terms  of  the  coefficient  matrices  of  the 
system,  as  in  the  following  theorems  [  3  ]  . 


Theorem  Z.  6:  A  fixed  eyetem  of  the  form  (2.  1)  ie  completely 


controllable  if  and  only  if  the  matrix  haa  rank  n  ,  where 

Q  =  [  B  :  AB  :  .  .  .  :  a"‘^b]  .  (2.  is) 

Theorem  2.  7:  A  fixed  system  of  the  form  (2.  1)  is  completely 
observable  if  and  only  if  the  matrix  has  rank  n  ,  where 

Qq  =  [c'  :  A'c  :  .  .  .  :(a""S'c']  .  (2.  i6) 

In  Chapter  3,  criteria  which  seem  to  be  natural 
generalizations  of  Theorems  2.  ^  and  2.  7  for  time-variable  systems 
will  be  derived. 


.  4^-  ‘  ■  ■*' " 
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CHAPTER  3 


CONTROLLABILITY,  OBSERVABILITY  AND  REDUCIBILITY  * 

3.  1  Introduction 

The  necessary  and  sufficient  conditions  for  controllability 
and  observability  summarized  in  Chapter  2  depend  explicitly  on  the 
system  solution  matrix  in  the  time-variable  case.  Since  this  matrix 
is  generally  not  available  except  as  a  numerically  tabulated  solution, 
it  has  been  difficult  to  characterize  the  structural  properties  of  a 
time-variable  system  in  terms  of  its  coefficient  matrix  description. 

In  order  to  circumvent  the  problem  of  system  solution 
we  propose  to  determine  the  extent  to  which  the  controllability  and 
observability  properties  can  be  characterized  in  terms  of  the  matrices 
A,  B  and  C  .  This  problem  has  been  studied  by  several  other 
workers  including  Stubberud  [  1] ,  and  Chang  [  2  ] .  Our  independently 
derived  results  will  be  described  below  and  compared  with  those  of 
the  above. 

Of  prime  importance  in  this  development  are  the  system 
"controllability"  and  "obser'^^ability"  matrices.  It  will  be  shown  that 
these  matrices,  which  do  not  require  knowledge  of  the  system  solution 
for  their  construction,  provide  significant  structural  information. 

Parts  of  this  chapter  have  appeared  in  two  papers  by  the  author  and 
H.  E.  Meadows  [  3,  4  ]  . 


25 


including;  (i)  a  sufficient  condition  for  complete  controllability  and 
observability,  and  (ii)  a  necessary  and  sufficient  condition  for  total 
controllability  and  observability.  These  conditions  are  established 
here  by  relating  the  controllability  and  observability  matrices  to  a 
new  test  for  linear  independence  of  certain  vector  functions.  This  test 
is  a  generalization  of  the  familiar  Wronskian  determinant  criterion 
for  scalar  functions. 

The  controllability  and  observability  matrices  will  also  be 
shown  to  be  useful  in  determining  whether  a  system  is  non- controllable 
or  non-observable  and  in  constructing  transformations  of  coordinates 
to  reduce  such  systems  to  zero- state  equivalents  of  lower  dynamic 
order. 

3.  2  The  Controllability  and  Observability  Matrices 

The  controllability  and  observability  matrices  of  a  time- 
variable  linear  system  of  the  form  (2.  1)  will  now  be  defined.  To 
insure  that  these  matrices  exist  and  are  well  behaved  for  all  time  it 
will  be  assumed  that  the  matrices  A,  B  and  C  are  continuously 
differentiable  n-1,  n  and  n  times,  respectively.  This  restriction 
can  be  considerably  weakened  for  many  argiunents  [  3  ] ,  but  at  the 
expense  of  greater  mathematical  complication. 

The  controllability  matrix  of  system  (2.  1)  is  defined  as 

Q  =  [Pn :  p .  :p^  ,)  o.i) 

c  U  .  1  .  .n-1 


■N 


where 


■*■  dt  ’ 


P,  =  B. 


(3.  2) 


The  observability  matrix  is  defined  in  a  dual  manner  as 


Oo  =  iSois,  ; 


;  ®n.l) 


(3.  3) 


=  A'Sk  "  5r®k  - 


Sfl  =  c. 


(3.4) 


Whenever  the  context  is  clear,  the  subscripts  'c'  and  'o'  will  be 


dropped. 


A  property  of  these  matrices  which  proves  to  be  quite 


useful  in  the  subsequent  development  concerns  their  behavior  under 
transformation  of  system  coordinates. 

A 

Property  3.  1 :  If  T  is  an  equivalence  transformation  and  if 

is  the  controllability  matrix  of  the  transformed  system,  then 


Q  =  TQ  . 
c  c 


(3.  5) 


then 


Equation  (3.  5)  follows  by  induction  since  if  P  = 


^  A  A  d  A 

^k+1  ^  ■‘^^k''’^^k 


=  (TA  +  T)T’^(TPj^)  +  (f  Pj^  +  Tlj^) 


=  T(-APj^  t  P^)  =  TP^^j 


1 


and  by  definition, 


P  =  B  =  TB  =  TP 

Q  ^  n 


Similarly,  if  is  the  observability  matrix  of  the  transformed  system 
then 


Qq  =  (T'l  Qj  . 


(3.6) 


To  indicate  the  role  played  by  the  controllability  matrix  as 
a  test  for  controllability  it  is  first  noted  that 


-^[x"^T)B(r)]  =  x"Sr)P  (T) 
dr^  ^ 


(3.  7) 


where  X(T)  is  any  fundamental  matrix  of  (2.  la).  Equation  (3.  7) 
follows  by  induction  since  P^  =  B  ,  and  if 


then 


1  1 

-^[X“  (t)B(t)]  =.x’  (T)Pj^(T), 
dT 


^j[x"  (r)B(T)]  = --[x'  (T)Pj^(r)] 

dT 

=  -X"  ^  (T)X(T)X'  ^  (T)  Pj^(T)+X‘  ^  (T)  Pj^(r) , 


• '  '-t 


L 


i 
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but 


X(T)X'Nt)  =  A(T)  . 


Therefore. 


,k+l  ^ 

— j^[x“  (t)B(t)]  =  x'  (T)[-A(T)Pj^(r)  +  p^(t)] 
dT 


=  x-\T)P^^,(r); 


It  follows  then  that 


[©(T) 


J3)'  '7  ' 


@ 


=  x"^r)Q  (r)  (3.8) 

c 


where  0(t  •  =  X  (t)B(t).  If  u  is  a  scalar  input  then  the  matrix  on  the 
left  side  of  (3.  8)  is  recognized  to  be  the  Wronskian  matrix  of  the 
rows  of  0(r)  .  It  is  well  known  [  5  J  that  the  rank  of  the  Wronskian 
matrix  may  be  used  to  test  the  linear  independence  of  scalar  functions. 
Since  X  ^  (T)  is  nonsingular  for  all  T  ,  the  rank  of  Q^('r)  is  equal  to 
that  of  the  Wronskian  matrix  for  all  T  .  For  single  input  systems, 
therefore,  the  controllability  matrix  yields  as  much  information  about 
the  degree  of  controllability  of  the  system  as  does  the  Wronskian 
matrix  of  the  rows  of  0(T)  .  In  the  following  section,  it  will  be  shown 
that  a  Wronskian  matrix  may  be  defined  for  vector  functions  having 
all  the  important  properties  of  the  scalar  Wronskian  matrix. 


/ 


■ 


V 


_ _ ,[ _ _  ....»>l■^'»*«':*'*n^i 


3.  3  The  Vector  Wronskian  Matrix 


Consider  the  set  of  r-dimensional  row  vector  functions 

9,{t/,  0,(t),  .  .  .  ,  0  (t)  ,  where  the  elements  of  each  0.(t/  together 
1  c  n  1 

with  their  first  n  derivatives  are  continuous  functions.  The  Wronskian 
matrix  of  such  a  set  of  functions  is  defined  as 


™1'92 . 9n>  =  C®0>i:  •  ■  • 


(3.9) 


where 


0Q  =  •  =0 


\  =  ■dr®k-i 


When  the  context  is  clear  W(t)  will  denote  the  Wronskian  matrix. 

The  following  theorem  provides  a  direct  generalization  of  the 
test  for  linear  independence  of  scalar  functions. 

Theorem  3,  1 :  If  W(t)  has  rank  n  for  some  t  in  an  interval  [t^.t^], 
then  the  functions  0,,  0_,  ....  0  are  linearly  independent  on  [t  ,  t  1- 
To  prove  the  theorem,  suppose  that  ti  e  0^  are  dependent 
on  >■  then  there  exists  a  constant  n-vector  \  0  such  that  for 

all  te[t^,  tj], 


X'0(t)  =  0 


.■i  ^  r  . 
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By  differentiattng  thi.  relation.hlp  n-1  times 


it  is  clear  that  for  all 


te  [tQ. tjl 


X'W(t)  =  0 


,hich  contradicts  W{t)  having 


rank  n  for  some  teCt^.t^]. 


The  followtag  theorem  provides  an  important  test  for  linear 
vector  functions.  Aithough  the  statement  of  the  theorem 
=  Vnovsn  result  tor  scalar  functions  [5]  , 


dependence  of  vector  functions.  Aitnoug..  u.. 

ig  .  direct  generalisation  of  a  hno»n  result  for  scalar  functions  t*]  , 

Che  proof  is  helieved  to  he  entirely  novel  in  that  it  utilises  properties 

corced  linear  systems  of  the  form  (i.  .)•  The  usual  proof  of  th.s 

_ rei  relies  on  the  relationship  the  Wronskian 


rcT  voties  on  the 

theorem  in  the  scalar  case  [S']  relies 

M  .  set  Of  functions  hears  to  the  solution  of  homogeneous  differential 
ecuations  -  no  such  analogous  relationship  holds  for  vector  functions. 
n.u.„rem  3.  Z'.  U  W(t)hasrank  q  less  than  n  for  all  te  [  t^,.  t,] 


tithe  wronskian  matrix  of  any  q  of  the  6,  has  rank  q  for  all 

fl  0  are  linearly  dependent  on 

then  the  functions  0^,  .  n 


t  1  and  may  be  expressed  as  a 


linear  combination  of  the  q 


independent  rows. 

Consider  the  system 


=  0(t)u 


(3.  10) 


r  -  t** 
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Let  Q  be  the  controllability  matrix  of  (3.  10);  then  it  is  clear  that  Q  =  W. 
Without  loss  of  generality,  it  may  be  assumed  that  6|  •  •  i  9  have 

a  Wronskian  matrix  of  rank  q  for  all  tj]  .  then 


where  has  q  rows  and  rank  q  for  all  te[t^,tj].  Since  Q  also 
has  rank  q  everywhere  on  the  same  interval,  the  rows  of  can  be 
written  as  a  linear  combination  of  the  rows  of  for  all  te[t^,t^]  : 
that  is 
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where  I  Is  the  q  "  order  Identity  matrix,  and  note  that 

q 


•‘V 

» 

.  0  - 


k  =  U,  1,  .  .  .  ,  n-1 


where 


^  dt  ^k-1  ’ 


'2  =0 


T  is  clearly  nonsingular  for  all  t  and  as  K  is  continuously  differentiable 


on  •  T  is  an  equivalence  transformation  on  the  interval.  If 


z  =  Tx  ,  therefore,  the  system 


z  =  Az  +  Bu 


y  =  Cz 


(3.13) 


where 


A  =  TT 


B  =  Tb 


C  =  T 


is  equivalent  to  (3.  10)  on  [  t^,  tj]  .  Since 


I  0 

q 


K  I 


A  = 


It  is  also  true,  however,  from  the  definition  of  P.  that 


A[P, 


“  “I  K; 


-i]  po  -  O,  ;■  •  •  :A„.,  -  n, 


-K  0  0  : .  .  . :  0 


Therefore, 


K  [fi. 


K  w(e,,  . e  )  =  o 


(3.  14) 


By  assumption,  W(9,,  82  , .  . .  8^)  has  rank  ,  for  all  teCt(,,tj]  so 
that  (3.  14)  implies 


K  =  0  , 


which  establishes  the  theorem. 


34 

The  following  corollary  to  Theorem  3.  2  provides  a  test  for 
determining  whether  a  set  of  functions  are  dependent  over  some  interval. 
Corollary  3.  1:  If  the  Wronskian  matrix  has  rank  less  than  n  for  all 
tc[to»ti]  then  there  exists  some  subinterval  of  over  which  the 

functions  •  ,6  are  linearly  dependent. 

16  n 

The  proof  of  this  corollary  by  induction  is  a  generalization 
of  that  for  the  scalar  case  given  by  Hurewicz  [5  ]  , 

For  n  =  1,  W(6j)  =  6j  so  that  if  the  rank  of  W(6j)  is  less 

than  one  for  all  then  6j  =  0  for  all  teJ^.tj]  and  the 

corollary  is  trivially  true. 

Suppose  that  the  corollary  is  valid  for  n  =  k-1  and 

W(8,,8, . 8.)  has  rank  less  than  k  for  all  t£[t  ,t,].  Either 

W(8,  ,8- . 8,  ,)  has  rank  k-1  for  some  te[t  ,  t.]  or  it  has  rank 

less  than  k-1  for  all  t€[t^,tj]  .  If  the  former  is  true,  then  by 
continuity  there  must  exist  some  subinterval  of  [t^,  t^]  on  which 
W(8, ,  6_,  .  .  .  ,  6  .)  has  rank  k-1  everywhere.  In  this  case,  Theorem 

3.  2  implies  that  0  ,,  0-,  ....  8,  are  dependent  on  the  subinterval.  If 

X  M  K 

W(0,,0,,...  0,  ,)  has  rank  less  than  k-1  for  all  te[t-,  t  ]  then  by 

assumption  0^,  0^,  •  ■  •  ,  6^  j  are  linearly  dependent  on  some  sub¬ 
interval  of  [  tp,  tj]  ,  which  completes  the  proof. 

It  is  now  possible  to  prove  the  major  result  for  the  vector 
Wronskian  matrix  which,  when  related  to  the  controllability  matrix, 
provides  a  necessary  and  sufficient  condition  for  total  controllability. 


* 


•>7 
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Theorem  3.  3;  A  necessary  and  sufficient  condition  for  the  row  vector 
functions  0^,  0^,  .  . .  ,  0^  to  be  linearly  independent  on  every  subinterval 
of  an  interval  [  t^,  t^]  is  that  their  Wronskian  matrix  not  have  rank 
less  than  n  on  any  subinterval  of  [1^,1^]. 

Sufficiency  follows  immediately  from  Theorem  3.  1  since 
any  subinterval  of  [  ^q>  niust  contain  points  where  W(t)  has  rank  n. 

Necessity  follows  from  Corollary  3.  1  since  if  W(t)  has 
rank  less  than  n  on  some  subinterval  of  [^Qjt^]  ,  there  must  be 
some  sub- subinterval  over  which  the  functions  are  dependent. 

If  the  elements  of  the  row  vectors  0,,  0_,  ....  0  are 
analytic  functions  of  t  the  Wronskian  test  for  linear  independence  can 
be  strengthened  considerably.  In  the  Appendix  properties  of  matrices 
of  analytic  functions  are  discussed,  and  the  following  resvilt  which 
follows  easily  from  Theorem  3.  2  is  established. 

Corollary  3.  2;  Let  0^,  0^,  .  .  .  ,  6^  be  row  vectors  of  functions 
analytic  on  [t-,t  ]  ;  then  0  ,  0  ,  . . . ,  6  are  linearly  independent 
on  [  t^,  tj]  if  and  only  if  W(t)  has  rank  n  for  some  t  c  [  • 

Furthermore,  if  W(t)  has  maximal  rank  q  less  than  n  ,  exactly 
q  of  the  0.  are  linearly  independent  on  the  interval. 

It  should  be  noted  that  the  rank  of  a  matrix  of  analytic  functions 
is  constant  save  possibly  at  a  finite  number  of  points  in  any  finite 
interval.  Thus  if  W(t)  has  rank  n  for  some  t  c  [  t^,  tj]  it  has 
rank  n  for  all  but  a  finite  number  of  points  in 
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[  t^] .  It  is  clear  then  that  if  a  set  of  analytic  row  vectors  are 

independent  on  any  subinterval  of  [  t^,  tj]  they  are  independent  on 
every  subinterval  of  • 


3.4  Criteria  for  Controllability  and  Observabilitv 


If  0(t)  =  X  (t)B(t)  where  X(t)  is  any  fundamental  matrix 


solution  of  (2.  1)  and  W  (t)  is  the  Wronskian  matrix  of  the  rows  of 

c 

0(t)  then  from  (3.  8) 


W^(t)  =  X“^t)Q  (t) 
c  c 


(3.  15) 


Immediate  application  of  Theorem  3. 1  gives  the  following  sufficient 
condition  for  complete  controllability. 

Theorem  3.  4;  (a)  System  (2.  1)  is  completely  controllable  on  the 

interval  [  t„,  t.l  if  Q  has  rank  n  for  some  te[t-,t,]. 

0  1c  0  1 

(b)  If  for  some  t  >  t.^,  Q  has  rank  n  then 

0  c 

system  (2.  1)  is  completely  controllable  at  t^  . 

To  simplify  notation,  if  an  analytic  matrix  has  rank  q  at  all 
but  a  finite  number  of  points  on  an  interval,  we  will  say  that  the 
matrix  "has  rank  q  "  on  the  interval. 
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This  result  was  also  prov' t.  directly  in  L3]  and  independently  by 
Stubberud  [1  ]  ,  and  Chang  [2],  and  is  the  time -variable  form  of 
Pontryagin's  "general  position  criteria"  [6].  It  was  also  shown  in  ’’3] 
that  Theorem  3.  4  holds  under  much  weaker  conditions  on  the  coefficient 
matrices. 

In  general  the  condition  of  Theorem  3.  4  is  not  necessary  for 
complete  controllability  as  demonstrated  by  the  following  example. 
Example  3.1:  Consider  the  second  order  single-input  system 


• 

*1 

■  0  0 

■^r 

g(t) 

z 

+ 

• 

*2-1 

.f(t)  0. 

-*2’ 

.  0  . 

u 


where 


0 

Iy(t-i) 
0 

y(t-5) 

0 

g(t)  =  f(t-2)  . 


t  s  0 

0  s  t  s  2 
2  :St  s  4 
4  ^  t  s  6 
6  s  t 


and 


exp( 


-1 


1  -  t‘ 


t  s  1 


y{t)  = 


t  il  . 


■I-., 


A  fundamental  matrix  for  this  system  is 


-r 


X(t)  = 


1  0 


[F(t)  1_ 

where  F(t)  =  J"  f(s)ds,  so  that 


®(t) 


(t)B(t)  = 


‘  1  0* 

g(t) 

g(t) 

_F(t)  1_ 

0 

_F(t)g(t)_ 

It  is  clear  here  from  Fig.  (3.  1)  that  g(t)  and  F(t)g(t)  are  linearly 
independent  on  the  interval  (0,8  ] .  However, 


g(t) 

g(t) 

g(t) 

|(t) 

_  0 

-f(t)g(t). 

_0 

0  . 

so  that  Q  has  rank  <2  for  all  t  e  [  0,  8  1 . 
c 

Observe  that  even  though  this  system  is  completely 
controllable  on  [  0,  8  ]  it  is  not  controllable  on  any  interval  [  0,  t  ] 
for  t  s  6.  That  is,  starting  at  time  t  =  0  there  is  a  waiting  period 
of  6  units  before  a  desired  final  state  can  be  reached. 

Although  a  necessary  and  sufficient  condition  is  not  possible 
for  complete  controllability  in  terms  of  the  matrix  such  a 

condition  can  be  given  for  the  more  useful  property  of  total  controll¬ 
ability.  This  condition,  which  follows  directly  from  Theorem  3.  3  and 
equation  (3.  15),  constitutes  a  major  resxilt  of  this  chapter.  It  may  be 
summarized  as  follows: 


I 


Theorem  3.  5:  System  (2.  1)  is  totally  controllable  on  the  interval 


[tQiti]  if  and  only  if  does  not  have  rank  less  than  n  on  any  sub¬ 
interval  of  [  t^,  tj]  . 


If  (2.  1)  is  a  fixed  system,  then  the  controllability  matrix 


takes  the  form 


Q  =  [ B  :  -AB  : 
c  . 


:(-a)''‘^b] 


and  is  c^uivalent  in  rank  to  the  familiar  form  of  the  fixed  controllability 


matrix  given  in  equation  (2. 15).  Theorem  3.  5  may  thus  be  viewed  as 


a  generalization  of  the  well  known  and  widely  useful  criterion  for 


controllability  of  fixed  systems  given  in  Chapter  2. 


It  should  be  quite  apparent  that  results  dual  to  those 


presented  for  control. ability  hold  for  observability.  To  formalize 


these  results  we  first  define  the  Wronskian  matrix  of  a  set  of  column 


vector  functions 


^(t)  =  [  4>  ^it)  >i>  ^{t) 


■0  (t)  ] 
•  n 


Wit)  = 


J 


(■ 


’/r 

y 


where 

*k'‘>  =  ^-^-1'“  -  =  *<*>  • 

If  ^(t)  =  C(t)X(t)  where  X(c)  is  a  fundamental  matric  solution  of 
(2.  1),  then  the  matrix 

WQ(t)  =  Q^X(t)  (3.16) 

plays  a  role  completely  dual  to  that  of  W^(t)  in  establishing  the 
following  theorems. 

Theorem  3.  6;  (a)  System  (2.  1)  is  completely  observable  on  the 

interval  [t^jt^]  if  has  rank  n  for  some  tC 

(b)  If  for  some  t  >  t^,  has  rank  n  then  the 
system  (2.  1)  is  completely  observable  at  t^  . 

Theorem  3.  7:  System  (2.  1)  is  totally  observable  on  the  interval  [  t^,  tj] 
if  and  only  if  does  not  have  rank  less  than  n  on  any  subinterval 
of  [  tg,  tj]  . 

If  the  coefficient  matrices  of  system  (2.  l)are  analytic 
♦ 

functions  of  time,  the  distinction  between  complete  and  total 
controllability  disappears,  as  will  now  be  shown.  First  note  that  if 
(2.  1)  is  an  analytic  system,  its  fundamental  matrices  must  also  be 
analytic  [  7  ] .  Thus,  it  follows  from  Corollary  3.  2  that  if  the  rows  of 

For  simplicity,  we  will  refer  to  (2.  1)  as  an  "analytic  system"  in  this 
case. 
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f 
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X  (t)  B(t)  are  linearly  independent  functions  on  some  interval  of  t  , 

they  are  independent  on  every  interval  of  t  .  That  is,  if  the  system 

is  completely  controllable  on  some  interval,  it  is  totally  controllable 

on  the  interval.  Since  the  converse  is  always  true,  complete  and 

total  controllability  are  equivalent  concepts  for  analytic  systems, 

and  controllability  is  independent  of  the  particular  choice  of  intervals 

and  initial  time  (as  is  true  for  fixed  systems).  As  a  special  case  of 

* 

Theorem  3.  5,  therefore,  we  have  the  following  corOi-.,ry. 

Corollary  3.  3:  An  analytic  system  of  the  form  (2.  l)is  completely 
controllable  if  and  only  if  has  rank  n  . 

Similarly  we  have: 

Corollary  3.  4:  An  analytic  system  of  the  form  (2.  1)  is  completely 
observable  if  and  only  if  has  rank  n  . 

3.  5  System  Reduction 

In  this  section  we  will  consider  the  problem  of  determining 
from  the  Q  matrix  whether  a  system  is  non- controllable  (i.  e. ,  not 
completely  controllable)  on  an  interval.  It  is  clear  from  Example  3.  1 
that  this  is  not  always  possible.  Coupled  with  this  problem  is  that  of 
"system  reduction."  It  is  well  known  [  8  ]  that  any  non- controllable 

This  result  was  also  established  by  Stubberud  [  1  ]  and  Chang  [  2  ]  . 


or  non-observable  system  can  be  reduced  to  a  zero -state  equivalent 
system  of  lower  dynamic  order.  When  possible  such  reduction  is 
usually  desirable  in  analysis  since  the  resulting  system  is  more  easily 
solved.  Moreover,  in  synthesis,  a  first  step  in  realizing  a  prescribed 
input-output  response  may  result  in  a  system  of  higher  order  than 
necessary.  A  least  order  realization  can  then  be  obtained  by  reducing 
the  initial  system.  This  type  of  synthesis  is  utilized  [  8  ]  for  fixed 
systems.  The  method  of  reduction  given  in{  8],  however,  is  based  on 
knowledge  of  the  transition  matrix  and  thus  cannot  readily  be  generalized 
to  time-variable  systems.  Reduction  of  time -variable  systems  is 
considered  in  [  9  ]  in  the  special  case  of  systems  realizable  by  a  single- 
input-single-output  differential  equations.  Presented  below  is  a  general 
approach  to  the  problem  of  system  reduction  based  on  properties  of  the 
controllability  and  observability  matrices.  The  following  will  be 
important  in  this  development. 

Property  3.  2:  If  has  rank  q  for  all  tc[  tj  ]  then  any  matrix  of 
the  form 


:Pj^]  ,  k>n-l 


has  rank  q  for  all  tc[tQ,tj]  . 

A  proof  of  this  property  is  given  in  the  appendix  of  [  3  ]  . 
Property  3.  3:  Let  be  the  controllability  matrix  of  (2.  1)  and  define 
the  related  matrix 


r 


.  - 

-♦ "  .  4; 


/ 


■'ih.  .1 


ifJinnu 


k 
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Qc  = 


(3.  17) 


Then  on  any  interval  where  has  rank  n  ,  the  system  A- matrix  can 
be  expressed  as 


•  ♦  t 

A  =  (Q  -  Q  )Q^ 
c  c  c 


Equation  (3.  18)  follows  from  the  observation  that 


(3.18) 


Q*  =  -AQ  +6 


c  c 


(3.19) 


Since  P.  =  B  ,  knowledge  of  consistent  matrices  Q  and  P  suffices 
0  "  n 

to  determine  (2.  1)  uniquely  on  any  interval  where  Q  bias  rank  n  . 

It  will  be  said  that  a  system  is  reducible  if  for  all  time,  its 
impulse  response  matrix  may  be  realized  by  z.  system  of  lower  order. 
Reduction  is  always  possible  if  a  system  is  not  completely  controllable 
and/or  not  completely  observable.  In  certain  degenerate  situations  a 
system  may  be  reducible  even  if  it  is  completely  controllable  and 
observable  [I'Jj  ,  but  such  cases  will  not  be  considered  here.  The 
reduction  problem  will  be  separated  into  two  parts; 

(i)  Construction  of  the  least  order  completely 
controllable  system  realizing  the  impulse 
response  matrix  of  (2.  1)  -  i.  e.  ,  that  part 
affected  by  the  input. 


‘a 


I 

1 


>|H 


\ 


_  _  jT 


(ii)  Construction  of  the  least  order  completely 
observable  system  realizing  the  impulse  response 
matrix  of  (2.  1)  -  i.  e.  ,  that  part  which  affects  the 
output. 

By  combining  (i)  and  {ii)  the  least  order  completely  controllable  and 
completely  observable  subsystem  may  be  found. 

Definition  3.  1:  System  (2.  1)  is  reducible  (from  the  input)  to  order 
q  £  n  and  no  lower  order  if  it  is  equivalent  to  a  system  of  the  form 


y  =  [Cj 


(3.20) 


and  the 


th 

q 


order  subsystem 


% 


B,u 


y 


(3.  21) 


is  completely  controllable. 

It  is  clear  that  the  impulse  response  matrix  of  (3.  21)  is 
equal  to  that  oi  (3.  20)  and  hence  to  that  of  (2.  1). 
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Definition  3.  2:  System  (2.  1)  is  reducible  (from  the  output)  to  order 
q  sen  and  to  no  lower  order  if  it  is  equivalent  to  a  system  of  the  form 


y  =  [Cj  0]z 


(3.  22) 


and  the 


th 

q 


order  subsystem 


"^ll'l 


+  BjU 


y 


=  c 


l"l 


(3.23) 


is  completely  observable. 

Note  that  in  this  case  the  reduced  system  is  zero- state  and 
zero  input  equivalent  (i.  e.  ,  "equivalent"  in  the  sense  of  [  1  1  ]  )  to  the 
original  system. 

Only  reduction  from  the  input  will  be  treated  below,  since 
reduction  from  the  output  may  be  effected  in  a  dual  manner. 

Conditions  under  which  a  system  is  reducible  (non-controllable) 
will  now  be  explored.  The  following  theorem  indicates  the  importance 
of  the  controllability  matrix  in  this  context. 


-u..- 


•v.  -M*'  ^  . 


Theorem  3.  8:  System  (2.  1)  is  reducible  to  a  q  order  totally 


controllable  system  if  and  only  if  an  equivalence  transformation  T(t) 


exists  such  that 


TQ  = 


(3.24) 


where  (Sj  has  q  rows  and  does  not  have  rank  <q  on  any  interval. 

A 

It  is  immediately  evident  that  if  Q  is  the  controllability 
matrix  of  a  system  of  the  form  (3.  20)  then 


Q  = 


(3.25) 


so  that  the  necessity  of  the  conditions  of  Theorem  3.  8  is  clear. 


To  prove  sufficiency,  it  is  first  noted  that  Property  3.  2 


implies 


where  Qj  has  q  rows,  so  that  from  (3.  19) 


„  *  A  d  • 

“^ll  "^12  °1  dT^l  ■  ^1 


^21  "^22  11° 


Since  Qj  does  not  have  rank  <q  on  any  interval,  and 
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^21^1  ^  ° 


it  can  be  concluded  that  A^j  =  d  .  Also,  since  S  is  formed  from  the  first  r 
columns  of  6  it  must  have  the  form 


B 


B. 


which  concludes  the  proof. 

A 

The  constraint  on  in  Theorem  3.8  can  not  in  general  be 
relaxed  to  include  reduction  to  systems  that  are  not  totally  controllable. 
It  is  possible,  as  shown  by  Example  3.  1,  that  a  system  may  possess 
a  controllability  matrix  of  the  form  (3.  25)  yet  not  be  reducible. 

The  conditions  under  which  Q  admits  a  transformation 
satisfying  Theorem  3.  8  will  now  be  examined. 

It  is  clear  that  if  such  a  transformation  T  exists,  then  its 
inverse  R  must  satisfy 


Q  =  RQ 


and  also  be  an  equivalence  transformation.  The  matrix  Q  may  be 


partitioned  as 
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where  is  a  q  x  (rq)  matrix  not  having  rank  <q  on  any  interval 

and  is  the  controllability  matrix  of  the  reduced  system  (3.  21).  Let 
R  be  partitioned  as  [  R  !  R  1  where  R  has  q  columns  and 
rank  q  everywhere.  Then  it  is  evident  that  a  necessary  condition  for 
T  to  exist  is  that  Q  not  have  rank  <  q  on  any  interval  and  that  the 

A 

first  rq  columns  of  Q  be  factorable  in  the  form  R  Q,,  .  If  such 

q  11 

factorization  is  possible  then  the  reducing  transformation  is  given  by 


T 


-1 


R 


(3.24) 


where  R  is  any  set  of  n-q  continuously  differentiable  columns 
n-q 

making  T  nonsingular. 

The  above  discussion  together  with  Theorem  3.  5  implies  the 
following  corollary  to  Theorem  3.8. 

Corollary  3.  2:  System  (2.  1)  is  reducible  to  a  totally  controllable 
system  of  order  q  ^  n  if  and  only  if  Q  does  not  have  rank  <q  on 
any  interval,  and  the  first  rq  columns  of  Q  can  be  factored  in  the 

A 

form  R  Q,,  ,  where  R  has  q  columns  and  rank  q  everywhere  and 
q  11  q 

A 

does  not  have  rank  <q  on  any  interval. 


N^. 


V 


mmmfnwitimmsttr-iriirptr.*'- 
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It  might  appear  that  the  conditions  of  Corollary  3.  2  are 
redundant  in  that  the  factorabiiity  of  Q  seems  to  be  implied  by  its 
rank.  That  this  is  not  the  case  is  demonstrated  by  the  following 
example. 

Example  3.  2:  Consider  the  second  order  system 


- 

■  • 

% 

fj(t) 

1 

where 


fj(t)  =  y{t) 
l^(t)  =  y{t  -  2) 

and  yft)  is  as  defined  in  Example  3.2.  For  this  system, 


Q 


t 


has  rank  1  almost  everywhere  on  the  interval  [-1,3],  as  does  the  first 
column  of  Q  .  The  system  is  completely  controllable,  however,  and 
is  thus  not  reducible  on  the  interval. 

A  simpler  sufficient  condition  for  reducibility  applicable  for 
a  somewhat  smaller  class  of  systems  is  provided  by  the  following 


corollary. 


5! 


Corollary  3.  3:  If  Q  has  rank  q  everywhere  and  q  columns  of  Q 
also  have  rank  q  everywhere  then  system  (2.  1)  is  reducible  to  a  q^*' 
order  totally  controllable  system. 

Any  set  of  columns  of  Q  having  rank  q  may  be  used  as 
in  (3.  24)^ avoiding  the  need  to  factor  Q  .  This  method  is  essentially  a 
generalization  of  Stubberud's  technique  (9]  for  single-input  systems. 
Cases  where  factorization  may  be  necessary  were  not  considered  in 
[9  J  ,  however,  nor  were  the  precise  conditions  under  which  reduction 
may  be  performed  made  clear. 

Observe  that  it  is  possible  for  the  Q  matrix  to  have  rank  q 
everywhere  without  any  q  columns  having  rank  q  everywhere,  as  in 
the  following  example. 

Example  3.  3:  Consider  the  system 


•  ■ 

^1 

sint 

Xj 

Sint 

L  2-1 

y  =  Xj  . 

The  controllability  matrix  of  this  system. 


sint 

cost 

sint 

cost 

has  rank  1  for  all  t  but  no  single  column  has  this  oroperty,  so  that 
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Corollary  3.  3  is  not  applicable.  It  is  clear,  however,  that  the  first 
column  of  Q  may  be  factored  as 


■r 

sint  , 

-1. 


so  that 


1  1 


1  U 


will  reduce  the  system  to 


z  j  -  (sint)u 


y  =  z, 


One  difficulty  involved  in  using  the  method  of  reduction 

associated  with  Corollaries  3.  2  and  3.  3  is  that  the  choice  of  the 

columns  R  may  not  be  obvious.  Theorem  3.  9  below  provides  a 
n-q 

method  of  reduction  in  terms  of  the  rows  of  Q  ,  which  has  the  virtue 
of  providing  a  simple  explicit  form  for  the  reducing  transformation  and 
the  reduced  system. 

Theorem  3.  9:  If  Q  has  rank  q  <  n  everywhere, and  q  rows  of  Q  , 
say  Qj  ,  have  this  property,  then  system  (2.  1)  can  be  reduced  by  the 
transformation 


✓  '  ■  1 


.  :,  r 
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T  = 


I  0 

q 


-K  I 


n-qj 


to  a  totally  controllable  system  of  order  q  where  the  rows  of  Q  have 
been  reordered  so  that 


Q  = 


Q 


1 


Q. 


and 


Q,  =  KQ, 


As  in  (3.  11),  K  may  be  expressed  as 


K  =  q^q;  , 


(}.  251 


so  that  T  is  an  equivalence  transformation.  It  is  clear  that 


TQ  = 


0 


which  proves  that  T  does  indeed  reduce  the  system. 

Although  (3.  25)  provides  an  explicit  form  for  the  proportionality 
matrix  K  ,  it  is  often  easier  to  find  K  by  alternate  methods.  A  method 
that  is  quite  useful  and  widely  applicable  is  presented  below. 


, 


I 


If  in  addition  to  the  conditions  of  Theorem  3.  9,  q  columns 


of  Qj,  which  will  be  denoted  as  ,  also  have  rank  q  everywhere 

then 


^  ^  ^21^11  ■ 


(3.26) 


where  Q  .  denotes  the  columns  of  Q  corresponding  to  the  columns 

b  1  M 

of  Qji  •  For  this  case,  simple  explicit  forms  for  the  matrices  of  the 
reduced  system  can  be  derived. 

Define  Qjj  to  be  that  submatrix  of  Q  which  corresponds 
to  Q u !  is, 


[°ul^  “u  • 

L^uJ 


Also,  let  Qj  be  the  first  q  rows  of  Q  . 


Then  since 


[:■] 


Q  ,  Q 

^11  ’  11 


it  follows  from  (3,  18)  that 


-  ..  .  ..  - ^ - 

f  ’  '  -r  \  -  a"  '  'V  .r  *S’. 

'Vz4  '  "i-  • 

■  V  1  ;S.'.  ‘.1  -  :  ^  *  ,  \  ^ 


1 


S>iJSC«iA9i^ 
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A 

A 


11 


-1 

11 


It  is  obvious  that 


A 

.  B, 


and 


C,  =  [C, :  c^] 


I 

q 

^21^11  . 


or 


=  C.  + 


C  Q  Q 
2^21^11 


-1 


(3.27) 


(3.28) 


(3.29) 


In  addition  to  giving  an  explicit  form  to  the  reduced  system  this 
reduction  method  has  many  computational  advantages.  Any  method  of 
system  reduction  involves  finding  a  transformation  of  coordinates.  The 
usual  reduction  procedure  is  to  find  the  matx'ices  of  the  reduced  system 
from  the  appropriate  submatrices  of  the  transformed  coefficient 
matrices  obtained  via  (2.  12).  This  approach,  which  has  been  used  by 
Kalman  [8  ]  for  fixed  systems  and  Stubberud  f9j  for  time -variable 
single-input  single-output  systems,  requires  the  inversion  and 
multiplication  of  n  x  n  matrices.  If  relations  (3.  26)  -  (3.  28)  are 
utilized,  however,  only  a  single  matrix  inversion  (of  a  q  xq  matrix) 
and  a  single  matrix  multiplication  (of  q  x  q  matrices)  are  required  to 

A  ^ 

find  Ajj  .  The  construction  of  is  also  simplified  and  no 

A 

computation  is  required  for  Bj  . 


\ 


\ 
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It  should  be  emphasized  that  this  explicit  method  of  system 
reduction  is  always  applicable  for  fixed  systems,  since  if  Q  has  rank 
q  it  is  trivially  true  that  some  submatrix  of  the  form  Qj^  also  has 
rank  q  . 


An  example  of  the  reduction  procedure  for  a  time-variable 
system  follows. 

Example  3.  4:  Consider  the  third  order  system 


' 

1 - 

o 

(N| 

o 

1 _ 

X 

1 

sint 

= 

-2  0  0 

^2 

+ 

cost 

.  ^3. 

cos2t  -sin2t  1 

.’‘3- 

sint 

The  controllability  matrix  of  this  system  is 


Q  = 


and  its  rank  is 


2 


sint 

-cost 

cost 

sint 

sint 

cost 

for  all 

t  .  Also 

-  sint 
-cost 

-  sint 

,  the  submatrix 


sint 

-cost 

cost 

sint 

/ 

/■  ■■  , 


v'-- 


/ 
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is  nonsingular  for  all  t  and  has  as  its  inverse 


j  rsint 

'll 

L-COSt 


sint  cost 
sint . 


For  this  choice  of  , 


j  =  [  s-int  cost] 


and 


Q 


11 


-cost 

sint 


-  sint 


-cost  J 


Therefore  it  is  seen  from  (3.26)  and  Theorem  3.9  that 


T  = 


0  n 


1  sinZt  1  J 


and  from  (3.  27),  (3.  28)  and  (3.  29)  that 


ro  2T  rsint -j  ^  ri  oi 

ol  J 


3.  6  Fixed  Systems 

Several  of  the  results  of  this  chapter  will  now  be  specialized 
and  applied  to  the  analysis  of  fixed  systems.  It  will  be  seen  both  here 


•'t, 

i 


t w ■  .If mutaeam 
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and  in  subsequent  chapters  that  the  more  general  approach 
necessitated  by  the  consideration  of  time -variable  systems  leads  to  a 
clarification  and  simplification  of  many  fixed  system  problems. 

The  reducibility  criteria  for  fixed  systems  simplify 
considerably,  as  evidenced  by  the  following  theorem 

Theorem  3.  10:  A  fixed  system  of  the  form  (2.  1)  is  reducible  from  the 
th 

input  to  a  q  -order  completely  controllable  system,  if  and  only  if 
has  rank  q  . 

The  proof  of  this  theorem  follows  directly  from  Theorem  3,8 

and  the  fact  that  if  Q  is  a  constant  matrix  of  rank  q  ,  there  always 

c 

exists  a  constant  non- singular  matrix  T  such  that 


TQ 

c 


0 


where  Q  ,  has  q  rows  and  rank  q  . 
c  1 

If  the  system  state-variables  are  so  ordered  that  the  first 
q  rows  of  Q  have  rank  q  ,  then 


I  0 

q 


-K  I 

n-qj 


where  K  is  given  by  equation  (3.  26),  and  the  reduced  system  by 
equations  (3.  28),  (3,  29),  and 
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Similarly,  we  have: 

Theorem  3.  11:  A  fixed  system  of  the  form  (2.  1)  is  reducible  from  the 
th 

output  to  a  q  -order  completely  observable  system  if  and  only  if 
has  rank  q  . 

An  important  advantage  of  the  criteria  for  reducibility  given 
by  Theorems  3.  10  and  3.  11,  as  well  as  the  reduction  procedure 
presented,  is  that  they  require  operations  solely  on  the  controllability 
and  observability  matrices.  Furthermore,  as  will  now  be  shown,  they 
lead  to  an  explicit  measure  of  the  degree  of  the  least  order  system 
obtainable  by  reduction  from  both  the  input  and  the  output  (i.  e.  ,  the 
least  order  system  realizing  the  system's  impulse  response  matrix). 
Theorem  3.  12:  Let  and  be  the  controllability  and  observability 
matrices  of  a  fixed  system  of  the  form  (2.  1),  then: 

(1)  The  system  is  irreducible  [  8  ]  (completely  controllable 

I 

and  completely  observable)  if  and  only  if  the  matrix  Q»Q  has  rank  n  . 

0  c 

(2)  If  Q„Q  has  rank  q  <  n  ,  then  (2.  1)  is  zero-state 

0  c 

equivalent  to  an  irreducible  system  of  order  q  , 

V  =  Av  +  Bu 

y  =  Cv  .  (3.  31) 

Furthermore, 


I 

I 


‘n- 

\ 


^  > 
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Q 'Q 

0  c 


CA 


ca’-* 


CA^-' 


[b:  ab 


(3.32) 


The  procf  of  part  (1)  is  easily  established.  If  has 

rank  n  ,  then  both  and  Q  must  have  rank  n  which  in  turn 

0  c 

Implies  the  systeni  is  irreducible.  If  Q 'Q  has  ’ank  <n  then  either 

0  c 


0  c 


has  rank  <n  ,  and  the  system  is  reducible  either  from  the 


input  or  the  output . 

Part  (2)  is  established  by  first  transforming  system  (2.  1)  to 
the  equivalent  form 


11 

0 


y  =  [Cj 


0 

A 

0 

0 


13 

22  "^23 


33 


43 


14 


24 


34 


44 


r~ 

B, 


2  '  ,  *'2 

'+  ' 


B. 

0 

0 


J  L  U 


u 


Cj  C^] z 


(3.  33) 


where  Zj  is  a  q -vector  and  the  subsystem 


^1  =  "'ll"!  "  ^1“ 


y  =  c,z, 
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is  completely  controllable  and  observable.  The  procedure  for 

constructing  this  equivalent  system  follows. 

Let  q  i  n  be  the  rank  of  Q  ,  then  by  Theorem  3.  10  there 
c  c 

exists  a  constant  nonsineular  transformation  T  such  that  if  w  =  T  x. 

®  c  c  ' 


w  = 


'^1 


W, 


11  12 


22J 


w 


1 


L'^ZJ 


B. 


1 


u 


y  =  [Cj  ]w  , 


and 


Q  =  T  Q 
c  c  c 


Q 


c  1 

0 


where  w,  and  Q  have  q  rows  and  Q  has  rank  cl.  If  Q 
1  cl  Cl  '^0 

is  partitioned  as 


A  A 


Q  '  =  r  Q  '  i  Q  ']  =  Q  ' T 
^0  ^^01  .02-'  ^0  c 


where  has  q  columns,  then  clearly 

01  c 


Q'Q  =  Q'Q  , 
0  c  01  cl 


(3.34) 


Since  Q  ,  has  rank  q  , 
cl  c 


A  A  A 


Q.,;  =  Q 'Q  Q  ,'(Q  ,Q  ; ) 

01  0  c  cl  cl  cl 


Q'Q  Q^, 
0  c  cl 


•I  n.*«»rtr»/;--.  »  •"'♦V  '  '  .  .  I 


which  implies  that  the  rank  of  is  ^  q  .  Equation  (3.  34),  however, 

A  A 

implies  that  has  rank  a  q  .  Therefore,  Q^j  must  have  rank  q  . 

It  follows  then  from  Theorem  3.  11  that  there  exists  a  constant 
nonsingular  matrix  T^  such  that  if 


•Tq  0 

W  , 

0  I 

L  n-q,,  J 


then  (2.  1)  is  equivalent  to  a  system  of  the  form  (3.  33).  Thus, 


^^01  :^02^  r  ^0 


0  “ 
I 

J 


=  tQo;:°;«o3:2o4^ 


where  Q^j  has  q  columns  and  rank  q 
If  is  partitioned  as 


0  IfQ  il 

cl 

1  0 

n-q„JL  J 


where  Q  ,  has  q  rows.then 
cl  * 


QnQ  =  Qni'Q  1 

0  c  01  cl 


.-I  ■' 


and  Q  must  have  rank  q  .  Making  the  identifications  v  =  Zj 
cl 

A  =  Ajj.  B  =  Bj.  C  =  Cj  and 


=  [  B  !  A  B !  .  . 
cl  •  • 

.  ;a”‘^b] 

• 

,  =  [  c  :  A  c  :  . 

01  . 

.  .  KA^'S'C'] 

completes  the  proof. 

Example  3.5:  Consider  the  system 


*1 

1 

.1 

1 

• 

0 

-  * 

1 

^2 

-1 

0 

0 

1 

*2 

+ 

2 

^3 

-1 

1 

-1 

1 

*3 

2 

.1 

2 

-3 

2 

2 

.  - 

y 

=  [  1 

1 

0 

.l]x 

. 

For  this  system 

'l  1  1  1* 

2  10-1 

Q  = 

c  2  10-3 


2  1 


0  -1 


1  1  0-1 


1  -3  4-1 


1  10-3 


1-3  4  -1_ 


Although  both  and  have  rank  Z  ,  their  product 


1111 


Q 'Q  = 

0  c 


1111 


1111 


1  1  1  1 


has  rank  1,  so  that  the  system  is  reducible  to  a  first  order  system.  To 
find  the  reduced  system  we  first  reduce  the  given  system  from  the 
input.  Note  that  the  submatrix 


■  [; :] 


of  Q  is  nonsingular.  With  this  choice  of  Q,,  , 
c  11 


■  c .:] 


so  that  from  (3.  30), 


•  V  "J 
<.3V-.  a, 


■  r,  .  ■  '-w  , 


‘\r: ' 
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i 


11 


1  0 

■  2  1 


and 


A 

3, 


Cj  =  [1  0] 


By  inspection,  therefore,  the  final  reduced  system  is 


V  =  V  +  u 


y  =  V 


r- 

i 


It  will  be  shown  in  Chapter  4,  that  the  form  of  the  product 
given  by  (3.  32)  leads  to  a  new  direct  method  for  constructing 
a  least  order  system  zero- state  equivalent  to  (2.  1).  This  procedure 
avoids  the  successive  transformation  required  by  Theorems  3.  10 
and  3.11. 

3.  7  Analytic  Systems 

Most  of  the  results  obtained  in  the  previous  section  for  fixed 
systems  generalize  directly  for  systems  with  analytic  coefficient 
matrices,  as  will  now  be  demonstrated. 

Theorem  3.  1  3:  An  analytic  system  of  the  form  (2.  1)  is  reducible 
from  the  input  to  a  order  completely  controllable  system,  if  and 
only  if  has  rank  q  . 


T 
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The  necessity  of  the  reducibility  criteria  is  obvious.  To 
prove  sufficiency  note  that  if  (and  thus  W^)  has  rank  q  , 
Corollary  3.2  implies  that  exactly  q  rows  of  X  St)B(t)  are 
linearly  independent.  Therefore,  there  must  exist  a  constant  non¬ 
singular  matrix  Tj  such  that 


T,  W 
1  c 


W 

c 

0 


where  has  q  rows  and  rank  q  .  Clearly,  T  =  T^X  is 

an  equivalence  transformation  and  from  (3.  15), 


W 

c 

• 

0 

Thus,  Theorem  3.  8  implies  that  the  system  is  reducible  to  a  q^^ 
order  completely  controllable  system. 

In  the  case  that  q  rows  of  Q  have  rank  q  for  all 
t  ,  the  explicit  form  of  the  reducing  transformation  given  in  Theorem  3.  9 
may  be  utilized.  If,  as  may  happen,  does  not  have  maximal 

rank  for  some  t  ,  the  factorization  procedure  may  be  necessary  to 
find  the  transformation. 

In  a  dual  manner  to  the  above  it  may  also  be  shown  that: 
Theorem  3.  14:  An  analytic  system  of  the  form  (2.  1)  is  reducible 
from  the  output  to  a  q^^-order  completely  observable  system  if  and 
only  if  has  rank  q  . 
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Since  the  proof  of  Theorem  3.  12  depended  solely  on 
Theorems  3.  10  and  3.  11  it  is  clear  that  a  completely  parallel  proof 
based  on  Theorems  3.  1  3  and  3.  14  may  be  given  for  the  following. 
Theorem  3.  1  5:  Let  and  be  the  controllability  and 
observability  matrices  of  an  analytic  system  of  the  form  (2.  1),  then: 

(1)  The  system  is  irreducible  if  and  only  if 

the  matrix  Q' Q  has  rank  n  . 

0  c 

(2)  If  QqQ^  has  rank  q  <  n  then  (2.  1)  is 
zero- state  equivalent  to  an  irreducible  system 
of  order  q  , 

V  =  Av  +  Bu 

y  =  Cv  .  (3.  34) 


Furthermore, 


Q'Q  - 

0  c 


n-  1 


] 


(3.  35) 


•»,  iP  ■  ‘i'- 


If 


i'^ 


}  .<• 


i 

I  ^ 
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In  Chapter  4,  it  will  be  ehown  that  an  extremely  efficient 
method  for  reducing  analytic  systems  follows  from  the  form  of 


given  by  (3.  35),  as  well  as  several  important  results  per- 


taining 


to  the  representation  of  such  systems. 


/  .  '  ' ''al-  .  ''i' i'  ja  ■  'Ji' 

fV. >■  .4*''- ■  i} 

iarfiV.^-  ••  ---a  e-  •  ^  I  I  \liMi  ■^*■1  1*  r**  I  ~~  fTTl  ^  *  • 
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CHAPTER  4 

SYSTEM  EQUIVALENCE  AND  CANONICAL  REPRESENTATION* 

4.  1  Introduction 

A  new  degree  of  controllability  and  observability,  stronger 
than  those  previously  defined,  will  be  introduced  in  this  chapter.  It 
will  be  shown  that  the  state  of  a  system  possessing  these  properties 
can  be  controlled  and  observed  "instantaneously.  "  The  motivation 
for  defining  these  properties  (to  be  called  uniform  controllability 
and  uniform  observability),  arises  from  the  role  they  play  in  the 
study  of  equivalent  systems.  There  are  many  problems  in  the  analysis, 
synthesis  and  control  of  linear  systems  that  can  be  solved  most 
effectively  by  transforming  a  given  system  to  an  appropriate  equivalent 
representation.  For  example,  when  simulating  a  time-variable 
system  on  an  analogue  computer,  it  is  desirable  to  minimize  the 
number  of  variable  components  (e.  g. ,  multipliers)  required.  There 
are  several  well  known  canonical  structures  for  single-input  single¬ 
output  systems  which  require  no  more  than  Zn  multipliers  [  1,  Z,  3]  . 

It  is  therefore  desirable  to  know  the  conditions  under  which  a  system 
has  an  equivalent  canonical  representation  and  to  be  able  to  construct 
such  representations  when  they  exist. 

For  fixed  systems,  conditions  for  a  system  to  be  equivalent 
to  several  particular  canonical  forms  are  well  known  [  1  ] ,  but  the 

Parts  of  this  chapter  have  appeared  in  two  papers  by  the  author  and 
H.  E.  Meadows  [  18,  19]. 
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techniques  for  transforming  to  these  forms  are  often  complex  or  ad-hoc 
methods  [  4,  5,  6] .  There  are  few  previous  results  for  time- variable 
systems  in  this  area. 

We  will  present  below  a  general  approach  to  the  problem  of 
system  equivalence  for  both  fixed  and  time-variable  systems.  For  the 
class  of  uniformly  controllable  and  observable  systems  (which  includes 
fixed  completely  controllable  and  observable  systems),  a  necessary 
and  sufficient  condition  for  equivalence  of  two  systems  will  be  given 
which  requires  no  knowledge  of  their  solutions.  In  addition,  an  explicit 
form  for  the  transformation  between  equivalent  systems  is  found.  New 
insight  into  the  problem  of  system  reduction  is  gained  from  this  study, 
and  a  method  for  directly  finding  the  class  of  minimal  order  systems 
zero- state  equivalent  to  a  prescribed  system  will  be  given.  An 
interesting  and  potentially  useful  byproduct  is  a  necessary  and  sufficient 
condition  for  an  analytic  system  to  be  zero- state  time -invariant  and 
a  method  for  finding  an  equivalent  fixed  system  if  such  is  the  case. 

Two  important  canonical  systems  will  be  studied  in  detail- 
the  input-output  differential  equation  form,  and  the  "phase- variable" 
canonical  form.  It  will  be  shown  that  uniform  observability  and 
uniform  controllability,  respectively,  are  necessary  and  sufficient 
for  the  existence  of  equivalent  systems  in  these  forms.  Explicit 
methods  for  constructing  canonical  equivalents  will  also  be  given. 


4.  Z  Uniform  Controllability  and  Uniform  Observability 


In  the  previous  chapter  conditions  for  complete  and  total 

controllability  were  given  in  terms  of  the  matrix  .  As  shown,  it 

is  not  necessary  for  this  matrix  to  have  maximal  rank  at  all  points 

to  insure  either  of  these  types  of  controllability.  A  point  to  be 

emphasized  in  the  present  chapter  is  that  for  many  problems  of  system 

equivalence  the  controllability  matrix  must  have  rank  n  for  all  t  . 

'Thus  we  define  the  following  new  degree  of  controllability. 

Definition  4.  1 :  System  (2.  1)  is  said  to  be  uniformly  controllable  on  an 

interval  [t^,  t,]  if  the  matrix  Q  has  rank  n  for  all  te[t.,,t,  ]. 

0  1  c  U  1 

A  new  degree  of  observability  is  defined  similarly: 

Definition  4.  2:  System  (2.  1)  is  said  to  be  uniformly  observable  on  an 
interval  [t^.t^]  if  the  matrix  has  rank  n  for  all  teLt^.tj], 

An  interesting  interpretation  of  uniform  controllability  can 
be  made  which  shows  how  this  criterion  relates  to  the  more  familiar 
types  of  controllability  which  arise  in  optimum  control  problems.  If  a 
system  is  totally  controllable,  then  by  definition  the  state  uf  the  system 
may  be  transferred  to  any  desired  value  in  an  arbitrarily  short  interval 
of  time  b/  application  of  some  input.  It  will  now  be  shown  that  if  a 
system  is  uniformly  controllable,  it  is  even  possible  to  perform  the 
state  transition  instantaneously.  Furthermore,  an  explicit  input  in 
terms  of  the  controllability  matrix  will  be  given  which  effects  the 


transition. 


It  is  first  necessary  to  examine  the  response  of  a  time> 
variable  system  to  impulse  functions  and  their  derivatives.  In  this 
discussion  we  will  rely  heavily  on  the  development  of  Zadeh  and  Desoer 
[3]  .  Since  the  response  to  an  impulse  is  discontinuous  in  general,  it 
is  necessary  to  distinguish  the  value  of  the  state  prior  to  and  after  the 
application  of  an  impulse.  The  notation  t  will  denote  the  left  hand 
side  of  any  discontinuity. 

Suppose  the  state  of  system  (2.  1)  is  zero  at  time  t^  and 
u(t)  =  6(t  -  1q)Q!q  where  is  an  r- vector  of  arbitrary  constants,  then 


1 

■I. 


»  (t,  t)B(t  )6(t  -  . 


‘>*0  • 


•  1 


By  the  sifting  property  of  the  impulse  function 


x(t)  =  $  (t, 


x(t)  =  $  [t.  tQ)B(tQ)aj^  , 


Thus,  it  may  be  said  that  the  state  at  time  t^  has  changed  instantaneously 


from  zero  to 


x(to) 


(4.1) 


,1*  $  ^ 


i 
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By  generalizing  the  argument  of  Zadeh  and  Desoer  ([3]  page  496)  to 

time-variable  systems  and  utilizing  (3.  7)  it  can  be  seen  that  if  u  (t)  = 

(k) 

6  (t  -  is  applied  to  system  (2.  1)  the  state  will  "jump"  to  the 

value 


*<‘o' 


^"o>“ 


k  ’ 


(4.  2) 


where  P,  (t)  is  defined  by  the  recursion  formula  (3.  2). 
k 

With  this  development,  it  can  now  be  shown  that  if  system 
(2.  1)  is  uniformly  controllable,  its  state  can  be  changed  instantaneously 
to  any  desired  value  with  an  input  of  the  form 

n- 1 

u(t)  =  IT  6^^^t  -  t^)a  .  (4.  3) 

k=0 

Let  X,  be  the  desired  state  at  t^  and  x.  be  the  initial  state 
d  0  i 

at  t”  ;  then  it  is  clear  that  if  u(t)  is  given  by  (4.  3) 

n-1 

'*d  -  ■‘i'  =  Ji 

k=0 

If  the  nr  vector  a  is  defined  as 


n,.-. 
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then  (4.  4)  can  be  rewritten  as 

(x  -  X.)  =  Q  (t^)Q!.  (4.  5) 

d  1  c  0 

Equation  (4.  5)  represents  n  equations  in  rr  unknowns,  and  if 

has  rank  n  ,  a  solution  (non-unique)  for  ot  exists.  An  explicit  solution 

is  given  by 

“=  =‘i> 

=  Qc<‘o’'''d  •  ’‘i>' 

If  the  matrix  A  (t)  is  defined  as 

A(t)  =  [6(t)i 

r  .  r.  •  r 

then  an  input  which  changes  the  system  state  from  x^  to  >■  ^  at  t^ 
is 

u(t)  =  A(t  -  t  -  X  )  .  (4.  6) 

U  c  U  Q  1 

For  a  single-input  system  the  above  solution  is  unique  and 
(4.  6)  reduces  to 


i,ffir^-tr^¥i  '^;*r>'tfir-r^‘&rni^-  J 


\  I 

*  # 

i  Si 


^(t  -  t(,)Q^  (l(,)(x^  -  X.) 


where 


C6(t)  :6^^^t) : 


It  is  interesting  to  observe  that  it  may  be  possible  to  change 


le  3Cate  instantaneously  at  t^  even  if  (^t^)  does  not  have  rank  n 


The  reason  for  this  is  that  appending  further  terms  of  the  sequence  of 


P  's  to  Q  can  increase  its  rank  at  a  point  even  though  it  cannot 

iC  c 


increase  its  rank  over  an  interval  (Property  3.  E).  A  simple  example 


which  demonstrates  this  is  the  first  order  system 


X  =  t  u 


For  this  system  =  t  ,  and  has  rank  0  at  t  =  0  .  However,  Pj  =  1 


and  the  input 


u(t)  =  6'  '(t)(x^-x^) 


will  change  the  system  state  from  to  at  t  =  0  . 


A  dual  result  for  uniform  observability  may  also  be 


established. 


Let  u  =  0  for  convenience,  and  differentiate  the  output  of 


(2.  1)  n-1  times.  Then  it  is  clear  that 
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y<^^t) 


Q^mit)  . 


(4.  7) 


If  QQ(tQ)  has  rank  n  ,  then  (4.7)  can  be  solved  uniquely  for  x(tjj)  , 
with  the  solution  given  explicitly  as 


or 

x(t(,)  =  (Q^(t|j))^Y(t|j)  .  (4.8) 

That  is,  if  system  (2.  1)  is  uniformly  observable,  the  state  of  the 
system  at  any  time  may  be  determined  instantaneously  from  obser¬ 
vations  of  the  system  output.and  its  derivatives. 

4.  3  Equivalent  Systems 

Consider  now  the  class  of  all  uniformly  observable  and 
controllable  systems  of  order  n  .  This  class  will  be  denoted  by  . 
If  a  system  of  the  form  (2.  1)  is  a  member  of  ,  the  notation 


(A.  B.  C)  e  U 

Cl 


will  be  used. 

If  (A,  B,  C)  and  (A,  B,  C)  are  equivalent  systems  of  order 

n  ,  where 


A  =  TAT"^  +  TT'^ 


B  =  TB 


C  =  CT 


(4.9) 


(4.  10) 


(4.11) 


the  relationship  between  the  two  systems  will  be  represented 
symbolically  as 


(A.  B,  C) 


If  (4.  12)  holds,  then  clearly 


(A.  B.  C)  . 


(4.  12) 


(A.  B.  C) - ^(A,  B.  C) 


(4.13) 


Two  equivalent  systems  (A,  B,  C)  and  (A,  B,  C)  will  be 
said  to  belong  to  the  same  equivalence  class.  The  primary  goal  of 
this  section  is  the  characterization  of  equivalence  classes  of  uniformly 
observable  and  controllable  systems.  It  will  be  shown  that  for  each 
such  class  there  is  a  unique  invariant  matrix,  knowledge  of  which 
suffices  to  generate  all  members  of  the  class.  We  will  use  this 
property  to  derive  several  important  results  including: 


-K.  I'-’  .  - 


■1 


t 

!■/ 
p  I 

%  I 
! 
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(1)  A  method  for  directly  finding  the  class  of  all  systems 
of  minimal  order  which  are  zero-state  equivalent  to  a  given  system. 

(2)  A  necessary  and  sufficient  condition  for  an  analytic  system  to  be 
zero- state  time -invariant  (i.  e.  ,  to  have  a  stationary  impulse  response 
matrix),  and  a  method  for  finding  an  equivalent  fixed  system  when 

such  exists. 

Recall  that  if  and  are  the  controllability  and 

observability  nnatrices  of  (A,  B,  C)  and  Q  and  Q  are  the 

c  0 

corresponding  matrices  for  (A,  B,  C),the  relationship 


(A,  B,  C) - ^(A,  B,  C)  . 


implies  that 


Q  =  TQ 
c  c 


Q '  =  Q  'T 


-1 


(4.  14) 

(4.  15) 


Clearly,  if  (A,  B,  C)  and  (A,  B,  C)  are  equivalent, then  (A,  B,  C)  e 

if  and  only  if  (A,  B,  C)  e  U  .  Equations  (4.  14)  and  (4.  15)  also 

n 

directly  imply  the  following  property  of  systems  belonging  to  U  . 

n 


Property  4.  1 :  If  (A,  B,  C)  e  U  and  (A,  B,  C)  0  U 

n  n 


are 


equivalent  where 
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m 

I 


< 

1 


(A,  B.  C) - 2_^(a.  B.  C) 


then 


T  =  Q 

c  c 


(Q  Q^)' 


(4.  16) 


or  equivalently, 


Q,,'Q  =  Q'Q  . 

0  c  0  c 


(4.  17) 


The  explicit  forms  of  the  transformation  between  equivalent 
systems  given  by  (4.  16)  will  be  extremely  useful  in  the  construction  of 
canonical  representations  for  time-variable  systems.  As  a 
preliminary  demonstration  of  the  utility  of  (4.  16),  it  will  be  used  to  derive 
an  explicit  expression  for  the  impulse  response  of  a  fixed  single-input 
system  (A,  b,  C).  For  simplicity,  it  is  assumed  that  A  has  distinct 
eigenvalues  .  Under  this  condition  it  is  well  known  [  1]  that  if 

(A,  b,  C)  is  completely  controllable  it  is  equivalent  to  (A,  b,  C),  where 


,  b  = 


The  controllability  matrix  of  (A,  b,  C)  is  seen  by  inspection  to  be  the 


Vandermonde  matrix  of  A  , 


1  X  j  .  .  .  \  j 

1  Xj .  .  .  x“‘* 


V  = 


1  X  .  . 
n 


It  follows  from  Property  4.  1,  therefore,  that 


(A,  b.  C) — I-MA.  b,  C). 


where 


T  =  VQ 


and  Q  is  the  fixed  form  of  the  controllability  matrix  of  (A,  b,  C), 


Q  =  [  b :  Ab :  .  .  .  :  A”"*b] . 


The  impulse  response  matrix  of  (A,  b,  C)  is  thus  given  by 


H(t)  =  CQV"  E(t)  , 


(4.  18) 


where 


V, 


-ftrmrt-'V.  ♦  \ifi| 
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E(t)  = 


Xjt 


\  t 
n 


Since  there  are  many  efficient  methods  for  inverting  Vandermonde 

matrices  [  7.8]  ,  equation  (4.  18)  would  appear  to  be  a  computationally 

advantageous  way  of  calculating  the  impulse  response  of  a  fixed  system. 

A  result  which  follows  immediately  from  Property  4.  1  is  a 

theorem  stated  by  Kalman  [1]  and  recently  established  by  Youla  [9]  . 

Theorem  4.1;  Let  (A,  B,  C)  e  U  and  (A,  B,  C)  e  U  be  fixed 
~  '  ■  n  n 

equivalent  systems.  Then 


(A.  B,  C) - ^—►(A.  B, 


where  T  is  a  constant  matrix. 

To  prove  this  theorem,  merely  note  that  Q  and  Q  must 
be  fixed,  and  apply  Property  4.  1. 

Although  Property  4.  1  gives  a  method  for  finding  the  form 
of  the  transformation  between  equivalent  systems,  we  do  not  as  yet 
have  a  criterion  for  determining  whether  two  given  sys  terns  are 
equivalent.  Before  considering  this  problem,  it  is  necessary  to 
introduce  two  new  system  matrices. 


s 

y.  *)■.'  ■ 


V  ’  .-A-v  . 


Define, 
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Q  =  L  :  p 

c  0.1 


•  n 


where 


and  define 


%  =  [Soi®! 


•s„^ 

.  n 


where 


S,  =  A'S,  ,  +  S,  ,  ;  =  C  . 

k  k-1  k-1  0 


*  + 

It  is  clear  that  and  are  proper  submatrices  of  ,  and  that 


Qp  and  are  similarly  proper  submatrices  of 


We  will  refer 


to  and  ot  as  the  identification  matrices  of  (A.  B.  C).  The 
c  U 

choice  of  this  name  is  made  clear  by  the  following  property. 

Property  4.  2:  The  matrices  and  B  (or  and  C)  are 

sufficient  to  uniquely  specify  (A,  B,  C)  e  U  . 

n 

Property  4.  2  was  essentially  proved  in  Chapter  3  since 


Wi 


.  «  + 
A  =  Q  -  Q  )Q^ 
c  c  c 


=  t(Q^C!o)Q;]'. 


(4.  19) 
(4.  20) 


and 


B  =  P. 


C  =  Sq 


Let  F  be  an  n  X  n  submatrix  o£  having  rank  n  for 

c 


X  A  T*  be  the  corresponding  submatrix  of 

all  t  (if  such  exists)  and  let  be  me  c 


Q  ,  that  is 

c 


,n,  equation  (4.  19)  may  be  aimpllfted  to 


«  #  -  1 
A  =  (F  -  F^)F^ 


(4.  21) 


Similarly,  let  Fq 


be  an  n  x  n  submatrix  of  Qj  having 


,an.  n  lot  all  t  and  tb.  corresponding  submatri.  of  (Q„) 


(4.  22) 


I.  is  now  possible  to  establish  a  necessary  and  sufficient 
condition  for  two  uniformly  observable  and  controllable  systems  to  be 
equivalent.  This  condition,  which  requires  no  Knowledge  of  the  system 
responses,  does  not  appear  to  have  been  shown  previously,  even  for 


fixed  systems. 


,  ■  ■ 


w- 


r 
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Theorem  4.  2:  Let  (A,  B,  C)  and  (A,  B,  C)  belong  to  ,  then  the 
two  systems  are  equivalent  if  and  only  if 


vx\va\  -  (q;)'q^ 


(4.  23) 


That  (4.  23)  is  necessary  for  equivalence  follows  directly 
from  Property  (4.  1). 

To  prove  sufficiency,  observe  that  if  (4.  23)  holds  then 

certainly 


=  Q' 

0  c  0  c 


(4.  24) 


since  this  equation  is  simply  the  first  ma  rows  of  (4.  23).  By 


assumption,  and  have  rank  n  for  all  t  ,  so  that 


q'^  =  tq'^ 

c  c 


where 


T  =  (Qi)^QJ 


P 


and  T  has  rank  n  for  all  t  .  Since  Q  has  rank  n  for  all  t  it 

c 

must  also  be  true  that 


(Qq)'  =  (qJ)'t"^  . 


?  •  •* 


. . .  _ 


r- 


w.  “  I 


•I  .1,  'V‘  "  I 
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The  form  of  T  implies  that  it  is  continuously  differentiable  and  thus 

A  A  A 

an  equivalence  transformation.  Consequently,  there  exist  (A,  B,  C)e 
such  that 


(A.  B.  C). 


(A.  6.  t)  . 


A  A  ^  A  A  A 

Let  Q  and  be  the  identification  matrices  of  (A,  B,  C);  then  (4.  14) 
c  0 

implies  that 


0“^  =  TQ'^  = 
c  c  c 


Thus  by  Property  4.  2, 

(A,  B,  C)  =  (A,  B.  C)  , 

and  (A,  B,  C)  is  equivalent  to  (A,  B,  C),  which  completes  the  proof. 

It  should  be  observed  that  it  suffices  for  equivalence  that 

either 


Wo>'Qc  = 


or 


Q  'Q^  =  Q  'Q^ 
0  c  0  c 


.  ...  .  ,  ’•  -..v  ■ 


The  matrix  T=  (Q^)'Q^  is  thus  seen  to  be  an  invariant 

0  c 

quantity  which  essentially  characterizes  an  equivalence  class  of 

systems.  It  will  now  be  shown  that  knowledge  of  T  alone  (and  not 

its  factors)  serves  to  determine  its  equivalence  class. 

Let  T  be  given  for  some  equivalence  class  of  systems  and 

let  (A,  B,  C)  e  U  be  any  member  of  the  class.  All  submatrices  of 

n 

T  arc  invariant  under  transformation  of  coordinates  so  that  we  can 
identify  tlie  following  useful  submatrices: 


(1) 


r=  Q'Q 
0  c 


is  the  matrix  formed  from  the  first  mn  rows  and  the  first  rn  columns 
of  T| 


(2)  r  = 


is  the  matrix  formed  from  the  last  mn  rows  and  the  first  rn  columns 
of  T. 

Suppose  that  F  contains  an  n  x  n  submatrix  F  having  rank  n 
for  all  t  ,  then  it  is  clear  that  F  is  of  the  form 


F  =  F-F 
0  c 


(4.  25) 


where  F^  and  F^  are  n  x  n  submatrices  of  Qj  and  respectively. 


J'..- 


Let  F*  be  defined  as  the  corre 


spending  submatrix  of  T  .  that  is 


F*  =  F*F 

*  0  c 


(4.  26) 


where 


*-0  =  "  *■» 


Also,  let  F,  be  the  n^ttlx  composed  of  those  columns  of  the  first 

„  ro»sof  r  which  correspond  to  F  .  that  is 


(4.  27) 


F,  =  CF 
‘  1  c 


iposed  of  those  rows  of  the  first 


Similarly,  let  F,  be  the  matrix  composes 

r  columns  of  T  which  correspond  to  the  rows  of  F  .  that  is 


(4.  28) 


IT  =  F  B  • 

*  2  0 


,  ■  F  F*  F  and  F,  are  all  submatrices  of 

Observe  that  the  matrices  F.  F  ,  *  j  2 

T  and  can  be  determined  ^dthout  knowing  the  matrices  of  any  particular 

j  ith  T  can  be  determined  directly  in 

equivalence  class  associated  with  T 
,  ns  of  the  above  submatrices. 

4  3.  LetTbegiven  for  an  equivalence  class  of  systems  in  , 

Theorem  4.  3.  ipeiiucg 

V*  F  d  F  be  defined  as  above.  Then,  the  system 
and  let  F.  F  ,  Fj  and  F^ 


(A,  B,  C)  is  a  member  of  the  equivalence  class  where 

-  ♦  -1 
A  =  F  F 


C  =  FjF"^  .  (4.29) 


Let  (A,  B,  C)  be  an  arbitrary  member  of  the  equivalence 

class.  Since  F  has  rank  n  for  all  t  then  both  and  F  must 

0  c 

have  rank  n  for  all  t  .  Consider  the  transformation  of  coordinates 


(A,  B.  C) 


(A,  B.  C)  . 


Under  this  transformation, 


and 


F  =  F  F  =  I 
^0  ^0^0  n  ’ 


F  =  F.F  , 
c  0  c 


A#  ♦  - 1 

^0  =  ^o*'o 


so  that  from  (4.  22) 
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''  -1  -1  - 
C  =  F.F  F_  =  C 
1  c  0 


Thus, 


(A.  B.  C) 


•^(A.  B,  C) 


and  the  theorem  is  established. 

The  fact  that  an  equivalence  class  of  systems  can  be 

determined  from  its  T  matrix  will  now  be  utilized  to  derive  a  new  and 

efficient  method  of  system  reduction. 

Let  (A,  B,  C)  be  a  fixed  system  of  order  n  with  r  inputs 

and  m  outputs  and  define  T  to  be  the  product  of  its  identification 

matrices  Q.  and  Q  .  If  T  has  rank  q  <n  ,  let  T  be  the  submatrix 
U  c 

formed  from  the  first  m(q+l)  rows  and  the  first  r(q+l)  columns  of 
T.  It  follows  then  from  Theorem  3.  12  that  there  exists  a  q^^-order 

AAA 

completely  controllable  and  observable  system  (A,  B,  C),  zero- state 

equivalent  to  (A,  B,  C).  Furthermore,  if  and  Q  are  the 

U  c 

identification  matrices  of  such  a  system  then 


T  =  . 


(4.  30) 


It  is  clear  that  Theorem  4.  3  may  be  applied  to  T  to  directly 

AAA 

construct  a  system  equivalent  to  (A,  B,  C)  (andthus  .-rero-state 
equivalent  to  (A,  B,  C)).  If  F,  F  ,  Fj  and  F^  are  defined  as  in 


« t.  ‘'Lf  'j.  ,  - 


■—•hjr'-:'-.  'T, 


Theorem  4.3,  as  submatrices  of  T  ,  then  the  following  theorem  holds. 

Theorem  4.  4:  If  (A,  B,  C)  is  a  time -invariant  system  and  T  has 

th 

rank  q  £  n  then  (A,  B,  C)  is  zero-state  equivalent  to  the  q  order 
completely  controllable  and  observable  system  (A,  B,  C)  where 


_  _  4c —  1 

A  =  F  F 


B  =  F. 


C  =  FjF"^  . 


(4.31) 


There  are  several  advantages  of  this  method  of  reduction 
both  in  comparison  with  previously  published  techniques  [  1,  lO]  and 


with  that  presented  in  Chapter  3.  Most  importantly,  reduction  from 
the  input  and  the  output  is  performed  simultaneously  with  a  resultant 
saving  in  computation.  Furthermore,  it  yields  a  simple  explicit  form 
for  the  reduced  system. 

To  compare  this  method  with  that  of  Chapter  3  consider  the 
system  of  Example  3.  5. 

Recall  that  for  this  system 


1  1 


Q 'Q  = 

0  c 


1  1 


1  1 


1  1 


/  *■  ^ 


'  ( •HMM’nr «  u  » 


JA*. 

f 


► 


V.- 


D 


♦  ' 

1  «»k 

'  I 

t 

i 
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Clearly,  QqQ^,  thus  T  )  has  rank  1  ,  and 


F  =  F*  =  F,  =  =  1  . 


Therefore, 


A  =  B  =  C  =  1 


Notice,  that  no  matrix  inversions  are  required  by  this  method  whereas 
by  the  previous  method  it  was  necessary  to  invert  a  2x2  matrix  to 
find  the  reduced  system.  It  will  generally  be  true  that  fewer  and  lower 
order  matrix  inversions  will  re  required  by  the  present  reduction 
procedure. 

An  obvious  generalization  of  Theorem  4.  4  to  analytic 
time- variable  systems  is: 

Theorem  4.  5:  If  (A,  B,  C)  is  an  analytic  system  and  T  has  rank 
q  s  n  for  all  t  then  (A,  B,  C)  is  zero-state  equivalent  to  the  qth 
order  uniformly  controllable  and  observable  system  (A,  B,  C)  where 


-  - 5l* —  -  1 

A  -  F  F 


B  =  F. 


C  -  FjF"\ 


(4.  32) 


. Wfc.Yt  WJ'i 
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Fj  and  F^  are  defined  as  in  Theorem  4.  3  as  submatrices 

It  should  be  noted  that  the  restriction  to  analytic  systems 
in  Theorem  4.  5  is  a  necessary  one  as  the  following  example  demon¬ 
strates. 

Example  4.  1 :  Consider  the  first  order  system 

X  =  g(t)u 
y  =  f(t)x 

where 

g(t)  =  y(t.i) 
f(t)  =  y(t-3) 

and  y(t)  is  as  defined  in  Example  3.  1.  It  is  clear  that  T  =  0  for 
t  .  However,  the  system  is  not  zero- state  equivalent  to  a  zero- 
order  system  since  its  impulse  response 

h(t,  T)  =  y(t-3)y(T-l)  ,  t  ^  T 

is  non  zero  for  T  c  [  0,  2]  and  t  e  [  2,  4] . 

An  interesting  and  potentially  useful  byproduct  of  Theorems 
4.  3  and  4.  5  is  a  test  to  determine  if  a  system  is  zero- state  time- 


and  F  , 
of  T  . 


■i 


iV 


i, 

(. 


i? 


i 


/ 


I 


I 

« 

i 

V 


invariant  (i.  e. ,  has  a  stationary  impulse  response  matrix),  and  a 
method  for  finding  a  fixed  zero- state  equivalent  for  such  a  system. 

This  effectively  provides  a  method  for  solving  a  class  of  time-variable 
systems  in  closed  form. 

Theorem  4.  6:  Let  (A,  B,  C)  be  an  n^^  order  time- variable  system. 
Then, 

(1)  (A,  B,  C)  is  equivalent  to  a  fixed  system  if  and  only 
if  T  is  a  constant  matrix  of  rank  n  . 

(2)  If  (A,  B,  C)  is  an  analytic  system,  it  is  zero- state 
equivalent  to  a  fixed  system  if  and  only  if  T  has  constant  rank  q 

s  n  and  the  submatrix  T  formed  from  the  first  m(q+l)  rows  and 

the  first  r(q+l)  columns  of  T  is  a  constant  matrix. 

The  proof  of  part  (1)  follows  from  Theorem  4.  3,  since 
He 

if  T  is  constant  then  F,  F  ,  and  F^  are  constant  matrices 
and  the  system  (A,  B,  C)  defined  by  (4.  29)  is  a  fixed  equivalent  of 
of  (A,  B,  C). 

Part  (2)  of  Theorem  4.  6  follows  similarly  from  Theorem 

4.  5. 

Example  4.  2:  Consider  the  time- variable  system 


/ 
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>;■> 

Also, 


F,  =  L  1  -1]  =  [  1  0] F 

V  ■ 

s,»  and 

f  r 


Therefore,  (A,  B,  C)  is  a  fixed  equivalent  of  the  above  system,  where 


4.  4  Canonical  Forms  for  Time-Variable  Systems 


The  transformational  properties  of  the  controllability  and 


,v  observability  matrices  will  now  be  utilized  to  determine  the  conditions 


,  I 


4 
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under  which  a  time-va.riable  system  has  equivalent  realizations  in 
several  canonical  forms.  Methods  of  transforming  to  canonical 
equivalents  will  also  be  given  when  they  exist.  Only  single-input- 
single-output  systems  will  be  considered  here,  and  to  avoid  the  chore 
of  counting  the  number  of  derivatives  required  in  various  arguments, 
it  will  be  assumed  that  all  system  coefficients  are  infinitely  differentiable. 

The  matrices  of  four  particular  canonical  structures  are 
given  below: 


_ 


0  a, 


0  a. 


.  0  a. 


:  .  b  = 


0  0  ...  1  a 


L  Cj  7!^  Cj  .  .  .cj 


(4.  34) 


0  .  .  .  0  a. 


0  ...  0  a^ 
c 


1  .  .  .  0  a 


’  ,b 


0  .  .  .  1  a 


[  0  0  0  .  .  .  Ij 


(4.  35) 


All  of  the  above  structures  have  the  advantages  of  requiring 
at  most  2n  time-variable  elements  for  their  realization.  Since  a 

general  single-input- single-output  system  of  the  form  (2.  1)  may  have 

?. 

as  many  as  2n  +  a  such  elements,  a  considerable  saving  of ’hardware 
may  result  if  a  system  can  be  realized  or  simulated  on  an  analogue 
computer  in  one  of  the  canonical  forms.  It  is  well  known  [  l]  that  any 
fixed  completely  controllable  system  has  equivalent  representations  in 


;.ri  '■ 


tn'f.'t 

'»\k 


(4.  34) 


(4.  35) 


All  of  the  above  structures  have  the  advantages  of  requiring 

at  most  2n  time-variable  elements  for  their  realization.  Since  a 

general  single-input- single-output  system  of  the  form  (2.  1)  may  have 

2 

as  many  as  2n  +  n  such  elements,  a  considerable  saving  of  hardware 


may  result  if  a  system  can  be  realized  or  simulated  on  an  analogue 
computer  in  one  of  the  canonical  forms.  It  is  well  known  [  1  ]  that  any 
fixed  completely  controllable  system  has  equivalent  representations  in 


i 


7; 
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forms  (4.  33)  and  (4.  34)  and  that  any  fixed  completely  observable 
system  has  equivalent  representations  in  forms  (4.  32)  and  (4.  35). 

As  will  be  shown  below,  however,  this  is  not  the  case  for  time-variable 
systems.  Only  uniformly  controllable  or  observable  systems  possess 
canonical  equivalents. 

Forms  (4.  32)  and  (4.33)  are  of  particular  interest.  System 
(4.  32)  is  the  standard  state-variable  representation  of  an  input-output 
differential  equation  [3]  ,  and  (4.  33)  is  the  "phase-variable"  canonical 
form  utilized  in  many  control  system  problems  [11-14].  The  former 
system  will  also  be  seen  to  be  valuable  in  synthesizing  prescribed 
impulse  responses. 

In  the  following  two  sections,  the  necessary  and  sufficient 
conditions  for  a  system  to  be  realizable  in  forms  (4.  32)  and  (4.  33)  will 
be  derived  and  methods  of  constructing  the  transformations  of 
coordinates  required  will  be  given.  As  forms  (4.  34)  and  (4.  35)  are 
dual  to  (4.  32)  and  (4.  33)  respectively,  it  is  not  necessary  to  treat  them 
separately. 

4.  5  The  Input-Output  Differential  Equation 

In  this  section,  the  conditions  under  which  system  (2.  1)  has 
an  equivalent  representation  of  the  form  (4.  32)  will  be  derived.  It  will 
first  be  shown,  however,  that  this  problem  is  the  same  as  finding  the 
conditions  under  which  (2.  1)  is  equivalent  to  a  differential  equation  of 
the  form 


■  h  'i  i 

I*',. 


,  V.. 


— irt  .f.j|-i - 
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n+1 

n 

E  v"'"  = 

E 

(4.  36) 

i=l 

U1  ' 

where  a  .  =  1  ,  and  the  coefficients  a.  and  b.  for  i  =  1,  2,  .  .  .  ,n 
n+i  1  i 

are  continuous  functions. 

The  equivalence  of  (4.  32)  and  (4.  36)  is  well  known  [  3]  and 
their  coefficients  are  related  by  the  equations 


a. 

1 


i  “  1 1  2|  •  .  .  )  n 


and 


b. 

1 


b  .  ,  + 
n-i+1 


n-i 


E 

j=i 


n-i-j+1 

E 


h-0 


h+i-1 

i-l 


)a 


i-j-h 


A 

b.  , 
J 


i  =  1,  2, .  .  .  ,  n-1 


b  = 
n 


A 

b 


1  • 


The  converse  problem,  that  of  t:xpressing  a  system  in  input- 
output  differential  equation  form  given  its  state-variable  representation 
has  been  considered  rather  intractable  by  several  authors  [  15,  16]  . 
Recently,  Weiss  and  Kalman  [  17]  have  presented  a  solution  which 
utilizes  the  impulse  response  of  the  given  system.  The  conditions  they 
derive  for  existence  of  the  input-output  differential  equation  are  some¬ 


what  imprecise,  however.  We  will  solve  this  problem  below  by  trans¬ 
forming  the  system  to  form  (3.  32)  which  in  light  of  the  above  discussion. 


1 


is  all  that  is  required.  Most  importantly,  this  solution  does  not  require 
knowledge  of  the  system  impulse  response. 

The  following  theorem  provides  a  necessary  and  sufficient 
condition  for  a  system  to  have  a  canonical  equivalent  of  the  form  (4.  3Z). 
Theorem  4.  7:  System  (2.  1)  (with  m  =  r  =  1)  is  equivalent  to  a  system 
of  the  form  (4.  32)  if  and  only  if  (2.  1)  is  uniformly  observable. 

Let  Qq  be  the  observability  matrix  of  (2.  1)  and  , 

A 

the  observability  matrix  of  (4.32).  By  inspection  Q.  =  I  ,  which 

0  n 

permits  the  following  compact  proof. 

Suppose,  (2.  1)  is  equivalent  to  (4.  32)  then  from  (4.  15) 


a  nonsingular  continuously  differentiable  matrix  T  exists  such  that 


Q  =  T'Q  =  T' 
'^0  0 


Therefore,  (2.  1)  must  be  uniformly  observable. 


Conversely,  if  (2.  1)  is  uniformly  observable  then  T  =  Qj 
qualifies  as  an  equivalence  transformation.  Under  this  transformation 
of  coordinates, 


A  d  '' 

A'  =  (-  ^Q  +  Q)Q 
'  dt  0  0  '  0 


Therefore, 


*  .  1  * 

O  -  O  O 
0  0  0  ■ 


where 


Also, 


c  = 


''°o' 


-1 


which  completes  the  proof. 

It  should  be  observed  that  the  above  proof  also  provide 
explicit  method  for  finding  form  (4.  3Z)  when  it  exists,  which  is 


summarized  below: 
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b  =  Q 'b  . 


Nota  that  contrary  to  the  results  of  [l?],  toUl  observability 

U  not  sufficient  for  the  existence  of  the  input  output  differential  equation. 
For  example,  the  system 


X  =  u 


y  =  (sint)x 

cannot  be  transformed  to  form  (4.  36). 

To  demonstrate  the  value  of  transforming  to  canonical  form 

lor  analogue  computer  simulation  or  analysis  »e  consider  the  following 

system; 

Example  4.  3 : 


cost  Sint 

2 

1  +  cos  t 

sint 

2 

-(1  +  sin  t 

-cost  Sint 

cost 

.{cost  +  2sint) 

sint  -  2cost 

-3 

[  cost 

-sint 

0 

The  observability  matrix  of  this  example  is 


r '*■  !  j  V; 


.'Wh'  ' 


"  4  !fA?l 


and  since  it  has  rank  3  for  all  t  ,  transformation  to  form  (4.  32)  is 
possible.  If  z  =  Q^x  ,  it  can  be  verified  that 


Thus  by  transforming  to  canonical  form,  we  obtain  a  fixed 
equivalent  in  this  case.  It  can  be  shown  in  general  via  Theorem  4.  4, 
that  if  a  uniformly  controllable  and  observable  system  possesses  a 
fixed  equivalent,  its  canonical  representations  will  be  fixed.  This 
provides  further  motivation  for  transforming  to  these  system  structux-es 
when  possible. 

4.  6  The  Phase-Variable  Canonical  Form 

The  problem  of  transforming  system  (2.  1)  to  an  equivalent 
phase-variable  canonical  form  is  a  considerably  more  difficult  problem 


i 

‘  k  * 


than  that  considered  in  the  previous  section.  The  reason  for  this  is 
that  the  controllability  matrix  of  (4.  33)  is  a  function  of  the  coefficients 
a^  .  In  the  fixed  case,  Kalman  [l  ]  has  shown  that  a  necessary  and 
sufficient  condition  for  pjoh  an  equivalent  to  exist  is  that  system  (2.  1) 
be  completely  controllable.  It  will  be  shown  here  that  under  somi  vhat 
more  restrictive  conditions  (uniform  controllability)  a  system  of  the 
form  (4.  33)  equivalent  to  system  (2.  1)  can  also  be  found  in  the  time- 
variable  case.  In  addition,  a  general  method  is  given  for  constructing 
both  the  equivalent  canonical  system  and  the  transforming  matrix  T  . 
The  specialization  of  this  method  to  fixed  systems  is  of  independent 
interest  as  it  proves  to  be  considerably  simpler  than  previously 
published  techniques  [4],  [s],  [e],  Johnson  and  Wonham  [4]  con¬ 
sidered  the  problem  of  finding  T  when  (2.  1)  is  fixed  and  A  has 
distinct  eigenvalue  s,  and  Mufti  [5]  extended  their  technique  to  the 
case  of  multiple  eigenvalues.  Further  refinement  of  the  methods  of 
[4]  ?nd  was  made  by  Chidambara  [bj  .  The  methods  of  constructing 
T  employed  in  W .  [5]  and  [6]  ,  however,  involve  finding  the  eigen¬ 
values  and  eigenvectors  of  A  ,  a  computational  task  equivalent  to 
solving  the  differential  equations  represented  by  system  (2.  1).  It  will 
be  seen  that  the  method  developed  here  is  independent  of  the  eigen¬ 
values  of  A  and  that  it  provides  simple  explicit  forms  for  the  matrices 
a  ,  T  and  T  ^  . 

Let  Q  ,  and  Q  be  the  controllability  matrices  of  (2.  1)  and 
(4.  33  respectively. 


I  -I  TiTTnilnr - - - - - -  (, 


Matrices  Q  and  p  will  be  examined  closely,  for  it  will 

n 

be  shown  that  they  serve  to  determine  the  transformation  from  (2.  1) 
to  (4.  33)  when  it  exists.  It  can  be  verified  by  direct  construction  that 


0  0  ...  0 


0  0 


(-1)  q 


n- 1 ,  n- 1 


^n,  n- 1 


Pn  = 


where 


0  -1 


‘^n-Z,  2  ‘^n-l.Z 


^  ^11  *  *  '  ^n-2,  1  '^n-l,  1 


(4.  37) 


^ik  ”  ^i-l,k-l  ^ '^i-l.k 


l<k<i<n 


-  (-l)'a 


i-2 

.  a  .q.  ,  .  ,  +  q.  ,  ,  ,  k=  l<i<n 

n-i+1  n-j  i-l,j+l  ^1-1,1 


=  (- 1)  a  ,  l<k  =  i<n  . 
n  "" 


(4.  38) 


From  the  form  of  Q  it  is  clear  that  any  system  of  the  form  (4.  33)  is 
uniformly  controllable.  A  more  informative  relation  for  q.^^  can  be 
easily  derived  from  (4.  38)  as 


108 


Vi+k 


i-k+1 

+  (-l)^  E  r  .q.  .  ... 

j=0 


+  2  (-D'^'^'q.  .  .  .^1  .  l<k<i<n 

i-j.k-j+1 


(4.  39) 


and 


.1  — 


q..  =  (-1)  a  ,  l<i<n  . 
ii  n 


It  follows  by  a  simple  induction  argument  that 


{terms  involving  only  the 

coefficients  a  , a  .  ,  , 

n  n-i+k+1 


! 


,  l£i<k<n 

(4.  40) 

For  notational  convenience,  the  bracketed  expression  in  (4.  40)  will 

be  represented  by  the  symbol  6  •  That  is,  any  function  that  can  be 

expressed  solely  in  terms  of  the  coefficients  a  , .  .  .  ,a  ,  will  be 

n  n-r+1 

replaced  by  the  symbol  wherever  no  other  information  about  the 
function  is  needed.  With  this  notation  (4.  39)  becomes 


*  "i-k  "  *i.ktl  ■  •  «o  =  0  • 

(4.4!) 

It  is  now  possible  to  prove 

Theorem  4.  8:  System  (2.  1)  is  equivalent  to  a  system  of  the  form  (4.  33) 
if  and  only  if(2.  1)  is  uniformly  controllable. 

The  necessity  of  the  controllability  condition  is  easily 
established,  since  if  (2.  1)  is  equivalent  to  (4.33)  where  z  =  Tx  ,  then 


' .  '  --  ’.iJ 


X'*  .“■••jL  ^ 


Q  =  TQ 


(4.  42) 


But  Q  and  T  have  rank  n  everywhere,  therefore  Q  must  have 
rank  n  everywhere,  which  implies  that  system  (2.  1)  is  uniformly 
controllable.  Conversely,  if  (2.  1)  is  uniformly  controllable  the  matrix 


T  =  QQ‘^ 


(4.  43) 


is  nonsingular  when  Q  is  the  controllability  matrix  of  any  system  of 
the  form  (4,  33).  Moreover  Property  4.  1  shows  that  (4.  43)  must  be 
the  form  of  the  transforming  matrix  if  it  exists.  Thus  to  prove  that 
the  uniform  -ontrollability  condition  is  sufficient  let  z  =  Tx  where  T 
is  given  by  (4.  43).  If  the  relationships 


A  =  (TA  +  T  )T 


-1 


(4.  44) 


and 


b  =  Tb 


(4.  45) 


can  be  satisfied  for  s  '.riii  un  ique  set  of  coefficients  a^  the  desired 
equivalence  wi’,<  ther.  be  eaiablished.  Note  that  (4.45)  is  valid  by 
inspection,  and  (4.44)  can  he;  rewritten  as 
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y  r'-'  ,  ' 


-  fi 


''y  - 


•S  •• 


-'  ■>  -Sali  -.i'wr  j 


'  f  •• 


i—  — d  ^ ^  “  1 

A  =  (QA  +  ^Q)Q 


(4.  46) 


in  which 


A  =  (Q  A 


■^A+  4-Q‘’)Q  =  -Q■^-AQ+  ^Q) 


-1  ♦ 
-Q  Q 


0  0  ...  0  g. 


.1  0  .  .  .  u  g. 


0  -1  ...  0  8: 


0  0  .  .  .-1  g. 


and  g  =  -Q-‘p„  .  where  i.  the  i-.h  element  of  the  vector  g  . 


From  (4.  46)  it  is  clear  that 


-aq+^Q  =  -QA 


■iPz-  •  •  Pn-l^nl  =lPjPj-  •  •  Pn-1<'°*>' 


Therefore,  it  will  suffice  to  satisfy 


Note  that  g  =  -'I  ,  as  defined  by  (4.  34). 


s£  i  C  'j 


i 


■i' 


A' 


«  IsiV'" 


n 


■Q  g  . 


(4.  47) 


It  follows  from  the  representation  of  the  elements  of  Q  and  p  given 

n 

in  (4.  37)  that  the  first  row  of  (4.  47)  is  a  =  g  and  that  remaining 

n  n 

rows  can  be  written  in  the  (symbolic)  form 


(“1)*^  +  ®  =  (-l)*^g  +  g  ,  1®,+*  •  .  g  6  ^  , 

'  r  n-r  ®r  ®r+l  1  *n  n-4  ’ 

-1  .  (4.48) 

Since  0^,  , . . ,  6^  ^  represent  terms  involving  only  the  coefficients 
^n  '  *  ’  ’  *^r+l’  determines  a^  uniquely  as  a  function  of 

the  known  g^,  •  •  • »  8^^  previously  calculated  a^  for  i>r  .  That 
is,  the  a^  can  be  found  recursively  beginning  with  “  8j^  • 
will  also  be  infinitely  differentiable  since  the  operations  involved  in 
(4.  48)  are  multiplication,  addition  and  differentiation  of  the  functions 
g^  which  are  infinitely  differentiable  due  to  the  assximptions  on 
system  (2.1).  The  matrix  T  found  by  substituting  the  recursively 
determined  a^  into  (4.  43)  is  therefore  nonsingular  and  continuously 
differentiable,  and  system  (2.  I)  is  equivalent  to  a  system  of  the  form 
(4.33). 

It  should  be  emphasized  that  the  proof  above  provides  a 
direct  method  for  finding  the  canonical  system  and  the  transforming 
matrix.  To  illustrate  the  method,  the  following  third  order  system 
will  be  considered. 


Ill 


f 

/ 


4  j  .^4  ■' 


'  S  • 


j>4r -.‘4^  V  . 


Example  4.4 :  Let  the  coefficient  matrices  of  system  (2.  1)  be 


1  1  0 


A  =  0  -I  -1  ,  b  =  0 

-e  ^  e  ^  -2  0 


The  matrices  Q,  p^,  and  g  are  then 


e  0  0  e  e 

^  n  n  -2t  -2t  -t 

Q  =  0  0  e  »P3=-  ®  »  8  =  e 

n  -2t  n  -3t  , 

0  e  0  e  1 


It  can  easily  be  established  from  (4.  48)  that  for  n  =  3 


^3  =*3-  *2“-e2  +  ^«3’  *1  =*1  -  *2  + *3 


so  for  this  example, 


0  0  1 


a  =  I  -e  I  and  Q  =  I  0  -1  1 


1  -1  1-e 


The  matrices  T  and  A  are  thus  given  by 


m  «<w«w  iwiiu  wiiiiiiaiiwiiiriilW 
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’o 

2t 

e 

0 

‘  0 

1 

o' 

T  = 

0 

2t 

e 

2t 

-e 

II 

\< 

0 

0 

1 

t 

,  e 

2t  t 
e  -e 

2t 

-e 

2.-'  -e-‘ 

1 

■ 

For  fixed  systems  the  transformation  method  outlined 
above  greatly  simplifies,  and  explicit  expressions  may  be  obtained 
for  both  the  coefficients  of  the  canonical  system  and  the  transforming 
matrix.  The  recurrence  relation  between  the  columns  of  the 
controllability  matrix  may  be  given  in  the  form 

Pk+l  =  ■^“k  >  Po  =  ^ 

and  from  the  fixed  equivalent  of  (4.  47)  it  is  found  that 

g  =  Q'^p^  =  Q'^A"b  (4.50) 

For  fixed  systems,  the  elements  of  both  g  and  a  are  the  coefficients 
of  the  characteristic  equation  of  A  so  that 

I  =  g  =  Q'^A"b,  (4.51) 


completely  specifying  the  canonical  system  (4.  33).  Once  the 
characteristic  equation  is  found,  the  transformation  T  is  given 

*  V 

As  are  a,  a,  and  a  . 


I: 

I 


"N 


f' 


_ 


'? 
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directly  and  explicitly  by  (4.43),  where 


0  .  .  .  0 


Q  = 


0-  ...  1 


C  o  c  -> 
n-3  n-2 


c,  .  .  .  c 


n-2  n-1 


(4.  52) 


and 


k-1 

=  Z  I 


'k  rr  '*n-iVi-l  ’ 
1=0 


k  =  l,2,...,n;  c.  =  1  . 


An  even  simpler  form  may  be  given  for  the  inverse  of  T 


-1 


(K  =  T  is  the  form  of  the  transformation  computed  in  references 
[  4j  -[  6] )  since 


-1 


-a. 


■a_  .  .  ,  -a 
3  n 


-a,  -a 


-a 


n 


1  .  .  .  0 


0  .  .  .  0 


as  may  be  easily  verified.  It  is  cleai  then  that 


rr 


I 
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K  =  T'*  =  QQ'^  (4.  53) 

requires  no  further  calculation  other  than  the  multiplication  of  Q  and 
Q  once  the  characteristic  equation  of  A  is  known. 

In  order  to  compare  the  computational  value  of  the  above 
procedure  with  existing  techniques  the  following  example  treated  by 
Mufti  [  5  ]  and  Chidambara  [  6  ]  will  be  considered. 

Example  4.  5 ;  Let  the  coefficient  matrices  of  system  (2.  1)  be 


1 

6 

-3  ■ 

"f 

A  = 

-1 

-1 

1 

,  b  = 

1 

.-2 

2 

0  _ 

1_ 

From  (4.  37) 


■  1 

4 

-2‘ 

‘10  ■ 

Q  = 

1 

-4 

-3 

’  P3  ’= 

-5 

1 

0 

-10_ 

-2 

The  coefficients  of  the  characteristic  equation  are  found  with  (4.  51)  as 


- 1 

o 

o 

1 

_ 1 

■  10  ■ 

'  -2 

7  -8  1 

-5 

= 

3 

1  4  -5 

-2 

0 

.jfc- 
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Thus 


so  that  from  (4.  53), 


"5f. 
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CHAPTER  5 

REALIZATION  OF  IMPULSE  RESPONSE  MATRICES 

5.  1  Introduction 

In  this  chapter  we  will  be  concerned  with  the  realization 
(or  synthesis)  of  prescribed  impulse  response  matrices  as  systems 
of  the  form  (2.  1).  Various  aspects  of  this  problem  have  been  treated 
by  Gilbert  [  1  ]  ,  Kalman  [  2,  3  ] ,  Weiss  and  Kalman  [  4  ]  ,  Youla  [  5  ]  , 
and  Desoer  [  6  ] ,  and  earlier  for  the  scalar  case  by  Darlington  [  7,  8  ] , 

Batkov  [  9  ]  and  Borskii  [  10  ]  .  We  shall  present  what  is  believed 
to  be  a  new  approach  based  on  the  properties  of  equivalent  systems 
derived  in  the  previous  chapter  and  on  a  method  of  Borskii  [  10  ] 
for  synthesizing  scalar  impulse  response  functions. 

Our  main  emphasis  will  be  on  the  realization  of  non¬ 
stationary  response  matrices  H(t,  T)  and  in  contrast  to  previous  work, 
we  will  not  assume  a  known  separated  form  for  H(t,  T)  (i.  e.  , 

H(t,T)  =  'Kt)3(T),  where  'l'(t)  and  ©(t)  are  finite  matrices).  The 
method  to  be  presented  here,  starts  with  a  non-separated  impulse 
matrix  and  yields  a  uniformly  controllable  and  uniformly  observable 
system,  when  such  exists. 

The  specialization  of  this  method  to  the  realization  of 
stationary  impulse  response  matrices  is  of  independent  interest.  It 
leads  to  an  easily  tested  necessary  and  sufficient  condition  for  a 
stationary  matrix  to  be  the  impulse  response  matrix  of  a  minimal  (completely 
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controllable  and  completely  observable)  fixed  system  of  specified  order, 
together  with  a  systematic  procedure  for  constructing  minimal  realizations. 

5.  2  Realization  of  Time-Variable  Impulse  Response  Matrices 

Let  H(t,  r  )  be  an  m  x  r  continuous  matrix  function  of  two 
variables.  It  will  be  said  that  H(t,  T )  is  realizable  if  there  exists  a 
system  (A,  B,  C)  for  which  H(t,  T)  lo  the  impulse  response  matrix. 

Any  system  having  H(t,  r  )  its  impulse  response  matrix  is  said  to 
realize  H(t,  T). 

It  has  been  shown  by  Kalman  [2]  that  a  necessary  and 
sufficient  condition  for  H(t,  r)  to  be  realizable  is  that  there  exist 
continuous  matrices  of  finite  dimension,  >Kt)  and  3(t)  ,  such  that 


k 


H(t,  t) 


'I'(t)3(T)  , 
0 


t  2T 

(5.1) 

t  <T 


The  proof  of  this  condition  is  quite  straightforward.  It  is 
clear  that  the  system  (0,  ©,  ^)  will  realize  (5.  1),  and  the  impulse 
response  matrix  of  any  system  of  the  form  (A,  B,  C)  is  given  as 


H(t,  T ) 


C(t)X(t)X 

0 


(T)B(T  ), 


tar 

t  <  T 


where  X(t)  is  a  fundamental  matrix  for  the  system. 

If  H(t,  T)  satisfies  (5.  1),  it  will  be  said  to  be  separable. 


u 

k 


The  major  problem  we  will  consider  in  this  chapter  is  the  realization 
of  response  matrices  which  are  not  given  in  separable  form,  or  even 
known  to  be  separable  a  priori.  There  are  several  known  procedures 


! 

} 

I 

I 

I 

I 


'<i' 
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for  realizing  fixed  systems  from  non- separated  response  matrices 
(usually  from  the  transfer  function  matrix  -  an  equivalent  specification 
to  a  non- separci ted  impulse  response  matrix)  [  1, 2,  3,  6]  ,  but  the 
general  time -variable  case  does  not  seem  to  have  been  treated 
previously.  We  will  present  a  procedure  which  starts  with  an  arbitrary 
matrix  function  of  two- variables  and  yields  when  such  exists,  a 
uniformly  observable  and  uniformly  controllable  realization.  The 
specialization  of  this  synthesis  procedure  to  fixed  systems  is  of 
independent  interest,  as  it  has  several  advantages  compared  with 
previously  published  techniques. 

Assuming  that  H(t,  T)  is  infinitely  differentiable  in  both 
of  its  arguments,  the  following  matrices  may  be  defined  for  all  i,  j  sO 
and  t  s  T  : 


3t^ 


H(t,T) 


.  .  .  Hp  ._j(t,r) 

H,o(t.T) 

HjjIt.T) 

•  •  • 

^..(t,T)  = 


^i-l,0^‘’^^  ^i-1,  1^^’’’^  •  •  ”i-l, 


th 


If  H(t,T)  is  realizable  by  an  n  order  system  (A,  B,  C)  then  it 
is  easily  verified  that 


/ 


I 


'I 


/ 
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H..(t,T)  =  S.'(t)f(t,T)P.(T) 
^  J  ^  J 


t 


where  $(t,T)  is  the  transition  matrix  of  (A,  B,  C)  and  S.  and  P. 

are  as  defined  in  Chapter  3.  Thus,  if  W„  and  W  are  the  Wronskian 

0  c 

matrices  of  (A,  B,  C)  then 


r  (t.r)  =  ^At)wlT) 

nn  0  c 


t  sT 


and 


r  (t.t)  =  Q;(t)Q  (t)  =  r(t,t) 

nn  u  c 


(5.2) 


Similarly, 


-rm-  <5.3) 

The  importance  of  Theorem  4.5  to  the  realization  problem 

is  now  clear.  Suppose  we  know  that  H(t,  T)  is  realizable  by  a  system 

of  order  n  (not  necessarily  a  least  order  realization).  The  T  matrix 
“th 

of  an  n  -order  realization  is  then  directly  obtainable  with  (5.  3)  so  that 
Theorem  4.  5  can  be  utilized  to  construct  a  uniformly  controllable  and 
observable  realization  if  such  exists. 

While  in  many  cases  it  is  possible  to  determine  n  by 
inspection  (this  is  true  for  all  stationary  impulse  response  matrices 


[2]  )  it  is  desirable  to  have  a  systematic  realization  procedure  which 

includes  the  determination  of  the  least  n  <  n  .  Observe  that  if  H(t,T) 

is  realizable  by  a  member  of  U  ,  there  exist  integers  tsn  and 

n 

k  s  ii  such  that  if  i  a  f,  and  j  then  1.  .(t,  T  )  has  rank  n  for  all 
t  i  T.  The  integers  and  k  are  the  least  such  that 


-1 


] 


and 


CP. 


have  rank  n  for  all  t  .  It  will  be  shown  below  that  under  certain 

conditions  the  converse  is  also  true.  That  is,  by  examining  the  ranks 

of  the  sequence  of  matrices  T  .(t,  t).  it  is  possible  to  determine  whether 

H(t,T)  is  realizable  as  a  member  of  U  .  Furthermore,  an  explicit 

n 

separation  and  realization  procedure  will  be  given. 

To  simplify  the  presentation,  we  will  first  consider  the 

case  where  H(t,  T  )  =  h(t,T),  a  scalar  function  of  two  variables. 

Theorem  5.  1:  If  for  some  finite,  non-negative  integer  n  the  matrices 

r  (t,  r)  and  r  ,  ,  (t,  T )  have  rank  n  for  all  t  and  r,  tar. 

nn  nt- 1 ,  nt-l 

then  h(t,  t)  is  realizable  as  an  n^^  order  uniformly  controllable  and 


uniformly  observable  system. 


ti-1 


h  .(t 

nj 


.T)  = 


i=0 


ai^j(t,  T)h(t,T)  ; 


j  =  0, 1, 


.  n 
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The  proof  of  this  theorem  together  with  a  canonical 

realization  of  h(t,  t  )  is  given  below. 

If  r  (t,  t)  and  i  .  1  “T)  have  rank  n  for  all  t^T  , 

nn  n+1, n+1 

then  certainly  the  last  row  of  F  .  ,  (t,  T  )  must  be  a  linear 

ht 1 , n+ 1 

combination  of  the  first  n  rows  for  all  t  and  T.  That  is, 


I' 


4' 

I 


for  some  set  of  coefficients  a.(t,  T  )  ,  i  =  1,  2,  . .  .  ,  n  .  As  will  now 
be  shown,  the  a^(t,  T)  are  unique,  continuously  differentiable  functions 
of  t  and  do  not  depend  on  T  .  Let 


and 


F(t,  T)  =  r  (t.r) 

nn 


F*(t,T)  =  ^F(t,T), 


A(t,  r)  = 


aj(t,  t)  a^(t,T) 


■v;:v 


.  0 

.  0 


(5.  4) 

(5.  5) 


.  a^(t,r) 


(5.  6) 


— ~~~r 

.f  V  'j*-  V  j.- k^'  V-, 


■A 


ir.v. 


V, 


T»  ^  ,  . 


-  ,  ,  •  I 


Since  F  (t,  T)  has  as  its  rows  the  last  n  rows  of 


r  .  1  (t,  r)  ,  it  is  clear  that 

n+ 1 ,  n  ’ 


F*(t,T)  =^F(t,r)  =  A(t,T)F(t,T) 
oi 


(5.7) 


Furthermore,  F(t, T)  has  rank  n  for  all  t  ^  T  so  that  equation  (5.7) 
possesses  the  unique  continuously  differentiable  solution 


A(t,  T)  =  F*(t,T)F'Nt,r). 


(5.8) 


To  establish  that  A(t,  T)  is  not  a  function  of  T  ,  differentiate  both 
sides  of  (5.  8)  with  respect  to  T  : 


^A(t,T)  =  [  ^  F*(t,T)]F’^t,T)  -  F*(t,T)F’\t,T)[  ^F(t,T)]  F"Nt,T), 

OT  cT  OT 

(5.  9) 

^  ^  jje 

Since  —  F(t,  T)  and  "^F  (t,  T)  are  formed  from  the  last  n  rows 

and  columns  of  F  ,,(t,  T)  and  F  ,,  (t,T)  respectively  it 

n,  n+1  ’  n+1 ,  n+1  ’  ’  '  ' 

must  also  be  true  that 


^  f\t)  -  A(t,r)  ^F(t,T). 

OT  OT 


(5.10) 


Therefore, 


•  w 


/ 


-^A(t,T)  =  A(t,r)[  ^F{t,T)]F'\t,r) 

-  A(t,  T)  F(t,  T)  F'  ^  (t,  T)[  ^  F(t,  T)]  F"  ^  (t,  T) 


=  0  . 


Thus,  with  a  slight  abuse  of  notation, 


A(t,T)  =  A(t) 


F*(t,T)F’^t,  T) 


(5.11) 


for  any  T  £  t  .  For  convenience,  let  T  =  t  in  (5,  11)  so  that  the 
defining  equation  for  A(t)  is 


A(t)  =  F*(t,t)F’Nt,t)  . 


(5.12) 


It  will  now  be  shown  that  F(t,  r)  can  be  separated  in  the 


form 


F(t,r)  =  F(t,  X)F“\t,  X)F(r,T)  . 


(5.13) 


for  any  arbitrary  parameter  X  . 

Observe  first  that  F  ^(t,  X)F(t,  T)  is  not  a  function  of  t  , 


since 


I' 


X)F(t,T)]  =  -F‘Nt,T)[-— F(t,  X)]F"Nt,  X)F(t,T) 

/ 

+  F(t,  X)~F(t,T) 

Ov 

=  -F(t,  X)A(t)F(t,  X)F'Nt,  X)F(t,T) 

+  F’^t,  X)A(t)F(t,r) 


/ 


One  realization  of  h(t,  T)  is  immediately  obtained  as 


z  =  F‘\t,  X)F(t,t)ej  u(t) 


y  =  e;F(t,  X)z 


(5.14) 


This  realization  is  without  feedback,  however,  which  for 
many  applications  may  be  undesirable.  Furthermore,  even  if  h(t,  r) 
is  stationary  (i.  e. ,  h(t,  T)  =  h(t  -  T)  ),  (5.  14)  will  be  a  time-variable 


system. 


A  more  useful  realization  is  obtained  by  observing  that 


for  an/  X,  F(t,  X)  is  a  fundamental  matrix  for 


^  =  A(t)x  , 


(5.15) 


where  A(t)  is  given  by  (5.  IZ).  Therefore,  the  transition  matrix  of 
(5.  15)  is  given  by 


0(t,T)  =  F(t,  X)F"  (T,X). 


The  impulse  response  h(t,  T)  can  thus  be  expressed  as 


(5.  16) 


h(t,T)  =  e;0(t,r)F(r,r)e, 


(5.17) 


and  Ib  clearly  realized  by  the  syctem  in  the  canonical  form  (4.  32) 


X  =  A(t)x  +  b(t)u 


y  =  cx 


(5.  8) 


where 


b(t)  =  F(t,t)e 


1 


A(t)  =  F*(t,t)F'Nt,t) 


(5.19) 


It  is  evident  from  (5.  2)  that  the  realizations  obtained  by 
this  method  are  uniformly  controllable  and  uniformly  observable, 
since  if  Q^(t)  and  QQ(t)  are  the  controllability  and  observability 
matrices  of  (5.  18)  then 


Oo'«  = 


and 


Q^(t)  =  F(t,t)  . 


An  important  property  of  the  realization  (5,  18)  is  that  if 
h(t,  T)  happens  to  be  stationary  the  system  will  of  necessity  be  time 
invariant.  This  is  obvious,  since  if  h(t,  T)  =  h(t  -  T)  then  it  must 


also  be  true  that 


F(t,  r)  -  F(t  -  T) 


and 


F*(t,  T)  =  F(t  -  T) 


so  that 


b  =  F(0)ej 


and 


A  =  F*(0)F"So)  , 


Example  5.  1  Suppose  we  are  given  the  scalar  function 


h(t,  r) 


e  sin(t-T), 


t  ^  T 


Observe  that 


r„(t,r)  = 


e  *ain(t-r) 


-t 

-e  cos(t-T) 


e  (cos(t-r).8in(t-r))  e  *(co8(t-T)+sin(t-T)) 


detr  (t,r)  =  e' 


t  s  T 


Also. 


r33(t,T)  = 


e  *sin(t-r) 


-e  ^co8(t-r) 


-t 

e  8in(t-r) 


e  (co8(t-r)-sin(t-r))  e  ^(cos(t-r)+8in(t-T))  e'^(co8(t-T)-sin(t-r)) 


-2e  co8(t-r) 


-2e  ^8in(t-T) 


•  2e  ^co8(t-r) 


and  clearly 


detr  (t,  T)  =  0, 


80  that  r  (t,  t)  also  has  rank  2  for  all  t  ^  T 


Let  F(t,  T)  =  r--(t,T)  ,  then 


1  32 


and 


F  (t,  t)  = 


-t  -t 

e  e 


-2e"^  0 


Applying  (5,19),  one  sees  that  h(t,  t)  is  realized  by  the  system 


It  should  be  clear  that  the  realization  procedure  above  is 
not  applicable  to  all  realizable  impulse  response  matrices.  If  a 
system  is  not  uniformly  controllable  and  observable,  or  reducible  to 
such  a  system,  its  impulse  response  matrix  will  not  satisfy  the 
conditions  of  Theorem  5.  1. 

Matrices  equivalent  to  F  (t)  and  F  ,  ,  ,  ,  (t)  were  first 

nn  nr 1 ,  nx 1 

used  by  Borskii  [10]  to  realize  single-input  single-output  systems. 
Borskii's  methods  and  proofs  are  somewhat  vague,  however,  and  lead 
to  several  false  conclusions  (such  as  a  statement  that  the  criteria  of 
Theorem  5.  1  are  also  necessary  for  realizability).  Furthermore, 
Borskii's  form  of  realization  is  an  input-output  differential  equation. 


which  result!  in  unnecessary  complication  in  finding  the  system 
coefficients. 


5.  3  Realization  of  Multi-Input  Multi-Output  Impulse  Response  Matrices 
We  will  now  consider  the  general  problem  of  realizing 
non- separated  matric  functions  of  two  variables.  The  major  result 
is  summarized  in  the  following  theorem. 

Theorem  5.  2:  Suppose  that  for  some  finite  positive  integers  1,  k 

and  n  (t  ^  n,  k  s  n)  ,  the  matrices  ^^+1  k+1^*’^^ 

have  rank  n  for  all  t  ^  T  ,  and  an  n  x  n  submatrix  of  P  ,  (t,  T)  also 

’  /,k 

has  rank  n  for  all  t  s  t  .  Then,  H(t,  r)  is  realizable  as  an 

i 

order  viniformly  observable  and  uniformly  controllable  system. 

A  constructive  proof  of  this  theorem  is  given  below. 

If  F(t,  T)  is  an  (n  x  n)  submatrix  of  T..  having 

^  JV 

rank  n  for  all  t  and  r  then  there  must  exist  matrices  Cj(t,  r) 
and  B^Ct,  T)  such  that 

H(t,T)  =  Cj(t,r)F(t,T)Bj(t,T)  .  (5.2->) 

Let  Fj(t,  r)  be  the  matrix  composed  of  those  columns  of  the  first 
m  rows  of  T  ,  (t,  t)  which  correspond  to  the  columns  of  F(t,  T)  . 

Then  it  is  clear  that 

Cj(t,T)  =  Fj(t,T)F'^t,r).  (5.21) 


Similarly,  let  F^(t,T)  be  the  matrix  composed  of  those  rows  of  the 
first  r  columns  of  r^j^(t,  T)  which  correspond  to  the  rows  of  F(t,  t)  . 
Then, 


Bj(t,r)  =  f"  (t,T)F2(t,r) 


Also,  let  F  (t,  T)  be  those  rows  and  columns  of 

corresponding  to  the  rows  and  columns  of  F  ,  (t,  t)  composing  F(t,  T) 

Vj  iC 

Then 


F*(t,T)  =  ^F(t,T)  =  A(t,T)F(t,r) 


As  in  the  scalar  case,  A(t,  r)  is  only  a  function  t  ,  and 


is  given  by 


A(t)  =  F*(t,t)F’St,  t) 


It  may  be  shown  similarly  that  Cj(t,  T)  is  not  a  function 


of  T  ,  nor  is  B(t,  T)  a  function  of  t  .  That  is, 


C,(t)  =  F,(t,t)F'  (t,t) 


B,(T)  =  F'  (T,T)F-(T,r) 


■W 

■I  Hill  ■  It* 


1^. 


Following  directly  from  (5.23),  as  in  the  single-input- 
single-output  case,  we  obtain  the  separated  form  of  F(t,  T), 


-1 


F(t,T)  =  F(t,  \)F  (r,X)F(r,T). 


(5.  26) 


Thus, 


-1 


H(t,T)  =  Cj(t)F(t,  X)F  (T,  X)F(T,r)Bj(r) 


(5.27) 


u 


t. 

I 

Si 

I" 


¥7: 

i'-} 


•'Ji. 


th 

and  is  realizable  by  the  n  order  system  (A,  B,  C),  where  A(t)  is 
given  by  (5.  23)  and 


B(t)  =  F2(t,t) 


-1 


C(t)  =  Fj(t,t)F  (t,t) 


(5.  28) 

(5.  29) 


Since  -f,  ^  n  and  k  s  n  ,  F  (t,  T)  must  have  rank  n  for  all  t  ^  T 

’  nn  ’ 


Therefore,  F  (t,  t)  has  rank  n  for  all  t  which  in  turn  implies 
’  nn  ’  ^ 


Q  (t)  and  Q  (t)  must  have  rank  n  for  all  t  .  The  system  (A,  B,  C) 

w  C 

A.  th 

*  is  thus  an  n  order  uniformly  observable  and  controllable 

realization  of  H(t,  r).  Furthermore,  as  in  the  scalar  case,  the  matrices 
A,  B  and  C  will  be  constant  if  H(t,  T)  =  H(t  -  r)  . 
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An  important  feature  of  this  realization  procedure  is  that 
the  systems  obtained  will  always  be  minimal.  A  system  is  said  to  be 


yrnm 


(1)  H(t,T)  =  >Kt)®(T)  , 


t  ^  T 


where  ^^(t)  and  ®(t)  have  n  columns  and  n  rows  respectively. 


(2)  H(t,T)  =  ^(t)®(T)  , 


t  i  T 


where  and  0(T)  have  n-1  columns  and  n-1  rows  respectively. 

If  r(t,  t)  is  calculated  from  (2),  it  clearly  has  rank  <n 
for  all  t  and  T  .  From  (1),  however, 


a  minimal  realization  of  an  impulse  response  matrix  H(t,  T)  if  there 
is  no  system  of  lower  order  which  also  realizes  H(t,  r).  It  is  clear 
from  our  earlier  discussion  of  reducibility  that  a  minimal 
realization  is  completely  controllable  and  observable  (globally  reduced 
[  5  ]).  The  converse  is  not  generally  true,  however, as  shown  by 
Desoer  [  6  ]  .  For  uniformly  (and  totally)  controllable  and  observable 
systems  the  converse  does  hold,as  will  now  be  demonstrated. 

Theorem  5.  3;  A  totally  controllable  and  totally  observable  system  is 
a  minimal  realization  of  its  impulse  response  matrix. 

Let  (A,  B,  C)  be  a  totally  controllable  and  totally  observable 
system  of  order  n  ,  and  suppose  that  a  system  (A,  B,  C)  of  order 
n-1  has  the  same  impulse  response  matrix  as  (A,  B,  C).  The 
impulse  response  matrix  can  then  be  written  in  two  forms: 


'ISBfRSltSjigSW!* 


r(t,T)  =  WjjWw^iT) 


t  2  T 


where  both  WQ(t)  and  W^(t)  do  not  have  rank  <n  on  any  interval. 

Consequently,  there  exist  values  of  T  and  t  with  t  5t  T  for  which 

r(t,  r)  has  rank  n  .  This  contradiction  establishes  the  theorem. 

Since  the  system  obtained  by  the  above  realization  procedure 

is  uniformly  controllable  and  observable,  Theorem  5.  3  implies  it  must 

be  a  minimal  realization  of  H(t,  T).  Furthermore,  it  is  unique  within 

an  equivalence  transformation  as  shown  by: 

* 

Theorem  5.  4:  If  (A,  B,  C)  C  i®  ^  realization  of  H(t,  T),  then  any 

other  minimal  realization  is  equivalent  to  (A,  B,  C). 

,  .. 

Let  (A,  B,  C)  beany  n  order  realization  of  H(t, T).  If 
Qq  ,  and  are  the  identification  matrices  of  (A,  B,  C) 

and  (A,  B,  C),  respectively,  then  clearly 


i  h 

I  r’ 


Thus,  it  follows  from  Theorem  4.  2  that  (A,  B,  C)  is  equivalent  to 
(A,  B,  C).  The  transformation  relating  the  two  systems  is  as  given  in 
Theorem  4.  2 


Theorem  5.  4  also  follows  from  Youla  ( [  5  ] ,  Corollary  1). 


5.  4  Synthesis  of  Fixed  Systems 


The  response  of  a  fixed  system  may  be  specified  by  its 
impulse  response  matrix  H(t  -  r)  or  its  transfer  function  matrix  G(s) 
(the  Laplace  transform  of  H(t)).  Realizability  of  H(t)  or  G(8)  is 
easily  determined  by  examining  the  form  of  each  element  of  the 
matrices  [  2  ]  .  For  example,  G(s)  is  realizable  if  and  only  if  all  of 
its  elements  are  the  ratios  of  polynomials  with  the  dominator  poly¬ 
nomial  of  higher  degree  than  the  numerator  polynomial  (i.  e.  ,  a  proper 
rational  matrix).  Finding  minimal  realizations  of  a  given  H(t)  or 
G(s)  is  far  from  a  trivial  problem,  however,  and  has  been  the  subject 
of  considerable  research  in  recent  years  [  1,  2,  3,  5,  6]  .  Several 
procedures  have  been  offered  for  constructing  minimal  realizations  and 
for  determining  their  order,  but  they  are  in  general  quite  complex  when 
G(s)  has  repeated  poles. 

We  will  present  below,  a  new  approach  to  the  realization 
problem  for  fixed  systems  by  specializing  the  results  of  the  previous 
section.  This  method  has  the  advantage  of  yielding  simultaneously, 
an  explicit  measure  of  the  degree  of  the  minimum  order  realizations 
of  H(t)  and  a  straightforward  synthesis  procedure  which  in  no  way 
depends  on  the  type  of  poles  possessed  by  G(s). 

Let  H(t)  be  an  m  x  r  matrix  function  of  t  .  Certainly, 
if  H(t)  is  to  be  the  impulse  response  matrix  of  a  fixed  system  it  must 
be  continuously  differentiable  an  arbitrary  number  of  times.  Thus 


we  may  define  for  all  i  ,  j 


H(t) 

H^^'(t)  .  .  . 

.  .  .  .  h0-'>(.) 

• 

H^^Nt)  .  .  . 

.  .  .  .  H^^t) 

H^‘Nt)  .  .  . 

where 


dt 


Save  for  the  sign  of  alternate  columns,  r..(t)  is  equivalent  to  r..(t,  0) 
as  defined  in  the  previous  section.  As  the  sign  changes  do  not  effect 
the  rank  of  any  of  the  matrices,  it  is  more  convenient  to  use  the 
present  definition. 

The  realizability  criteria  for  fixed  systems  may  now  be 
given.  The  proof  closely  follows  that  of  Theorem  5.  2. 

Theorem  5.  5:  H(t)  is  the  impulse  response  matrix  of  a  fixed  completely 

controllable  and  completely  observable  system  of  order  n  if  and  only 

if  there  exist  non-negative  integers  n  and  k  ^  n  such  that  the 

matrices  k+1^^^  have  rank  n  for  all  t^O  and  an 

n  X  n  submatrix  of  r  .  (t)  has  rank  n  for  all  t  ^  0  . 

ik 


wnm'Wup 
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Let  F(t)  be  an  n  x  n  submatrix  of  F  ,  (t)  having  rank  n 

9|e 

for  all  t  i  0  .  With  respect  to  this  matrix,  let  the  matrices  F  (t)  , 

Fj(t)  and  t)  be  defined  precisely  as  in  the  proof  of  Theorem  5.  2. 

Then, 


H(t)  =  Cj(t)F(t)Bj(t) 


(5.  30) 


where 


Cj(t) 


Fj  (t)F"^t) 


and 


Bj(t) 


=  F 


(5.31) 


(5.  32) 


Also, 


F*(t)  =  F^^^t)  =  A(t)F(t)  , 


(5.  33) 


since  both  F.,.  ,  , ,  (t)  and  F.  ,  (t)  have  rank  n.  Furthermore, 
/.+  l,k+l  I,,  k  ’ 

A(t)  is  not  a  function  of  t  since 


^A(t)  =  F^^^t)F’\t).  F^^^(t)F’Nt)F^^Nt)F"\t)  , 


:v 


and 

=  A(t)F^^Nt) 


imply  that 


=  0  . 


It  may  be  shown  ■ 


imilarlythat  C^U)  and  are  constant  matrices. 


Thus, 


F(t)  =  AF(t)  , 


(5.  34) 


and 


H(t)  =  CjF(t)Bj 


(5.  35) 


Since  F(t)  ha.  rank  n  for  all  t  ,  (5.  34)  implie.  that 


F(t)  =  e^*F(0) 


and  is  a 


fundamental  matrix  lor  the  .y.tem  of  n  differential  equation. 


X  =  Ax  . 


......a...  •mumuEdiMdi'u.  '  I-. jn, 

aiwa  Jiir 


.  -ffe-  - 


From  equations  (5.  31),  (5.  32)  e.nd  (5.  34) 


A  =  F(O)F'So) 

Bj  =  F'\o)F^{0) 

Cj  =  Fj(0)f‘No)  . 

Thus, 

- 1  At 

H(t)  =  Fj(0)F  (0)e  F^CO) 


(5.  36) 

(5.  37) 

(5.  38) 


(s.  39 


and  is  realized  by  the  n  order  system  (A,  B,  C)  where 


B  =  F^{0)  (5.40) 

and 

C  =  Fj(0)F‘\0)  .  (5.41) 

Furthermore,  if 

Q  =  [  B  :  AB  :  ...  :  a"”^b] 

0  •  •  • 


% 


v 


Vi' 


* 


-  Jr^  ».'  ’ 


v  ^  . 
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and 


Qq  =  [c  :a'c'  : 


.  :  (a""S'c'] 


then 


r 

nn 


(0) 


Q'Q 

0  c 


Since  P  (0)  has  rank  n  ,  (A,  B,  C)  is  a  completely  controllable  and 
nn 

observable  system. 

The  converse  of  Theroem  5.  5  is  easily  established.  If 
(A,  B,  C)  is  an  n^^  order  realization  of  H(t)  then  it  may  be  verified 
that 


r 

nn 


(t) 


At^ 

Q>  Q 


If  Q/i  and  Q  have  rank  n  ,  then  certainly  P  (t)  has  rank  n  for  all 
0  c  ’  ^  nn 

t  .  Similarly, 


has  rank  n  for  all  t  .  Also,  if  Q  and  Q '  have  rank  n  they  must 

c  U 

contain  n  x  n  submatrices,  and  ,  respectively,  having  rank  n. 


X 


Ir  rmamM 


li. _ 


Consequently, 
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F(t)  =  F-e  F 
U  c 


is  an  n  X  n  submatrix  of  T  (t)  having  rank  n  for  all  t  . 

nn 

It  follows  immediately  from  Theorem  5.  3  that  the  system 
(A,  B,  C)  derived  above  is  a  minimal  realization  of  H(t)  and  from 
Theorem  5.  4  it  is  assured  that  (A,  B,  C)  is  unique  within  an 
equivalence  transformation.  Furthermore,  Theorem  4.1  implies 
that  all  constant  minimal  realizations  of  H(t)  are  related  to  (A,  B,  C) 
by  a  constant  transformation. 

Example  5.  2;  Consider  the  matric  function 


H(t)  = 


-t 


e’^8+3t) 


2e  ^  e  ^(4+6t) 


-t 


The  determinant  of  H(t)  is  -12e  ,  so  that  Tj  j  (t)  =  H(t)  has  rank  2 


for  all  t  .  ^22^^^  found  to  be 


^22^’'^ 


-t 

e 

e‘^8+3t) 

-t 

-  e 

.e"^5+3t) 

2e'^ 

e  ^(4+6t) 

-2e"^ 

-e‘^-2+6t) 

-t 

-  e 

.e"^5+3t) 

-t 

e 

e’^2+3t) 

1 

(M 

-e“^-2+6t) 

2e'^ 

e’^(-8+6t) 

Denote  the  columns  of  ^2*  83  ^4  >  clear 


•  •  w  .t*.  ^  ir^  : 


,  V’’ 


that  g^  =  -gj  and  g^  =  3gj  “  82  »  *°  ^ 

aU  t  .  Let  F(t)  =  H(t)*  then 

Fj(t)  =  F^{t)  =  F(t) 


and 


F(t) 


=  H(t) 


Therefore, 


and 


F(0)  = 


F(0)  = 


-1 


-2 


=  Fj(0)  =  F^{0) 


•ii 

J 


Thus  H(t)  is  realized  by  the  second  order  completely  controllable 
and  observable  system  (A,  B,  C),  where 


-  =  -i 


2 

-4 


1 

6 


t 

t 


s' 

i' 


\ 


•/? 


B  = 


1  8 

2  4 


and 


C  = 


0 

1 


Although  the  above  procedure  provides  a  self-contained 
theory  for  minimal  realization  of  fixed  impulse  response  and  transfer 
function  matrices,  there  are  several  alternate  methods  which  might 
be  better  suited  for  synthesizing  high  order  systems  on  a  digital 
computer.  Among  these  are  Gilbert's  method  for  transfer  function 
matrices  with  simple  poles  [  1  ],  and  Oesoer's  extention  for  multiple 
poles  [  6  ]  .  The  advantage  of  these  methods  for  fixed  systems  is  that 
they  do  not  require  operations  on  time- variable  matrices.  We  will 
describe  another  such  computational  procedure  below  which  utilizes  the 
explicit  method  of  system  reduction  presented  in  Chapter  4. 

As  observed  in  the  discussion  prior  to  Theorem  5.  1,  one 
can  always  determine  by  inspection  the  order  (say  n  )  of  some 
realization  of  a  fixed  impulse  response,  or  transfer  function  matrix. 
Kalman  [  2  ]  uses  this  information  in  a  synthesis  procedure  which 
entails  first  constructing  a  realization  of  order  n  ,  and  then  reducing 
it  to  minimal  order. 


-  ‘i  — .  *1 

*  i 
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A  modification  of  this  approach,  which  will  now  be  outlined, 
greatly  reduces  the  computation  required  and  results  in  an  explicit 
form  for  a  minimal  realization  of  G(s)  .  This  method  is  based  on  the 
procedure  given  in  Chapter  4  for  obtaining  a  least  order  representation 
of  a  system  from  its  T  matrix. 

(1)  Determine  n  ,  an  upper  bound  on  the  order  of  the 
minimal  realizations,  by  inspection  of  G(s)  (one  such  method  for 
finding  n  is  given  by  Kalman  [  2  ] . 

(2)  Form  F _ (0).  It  is  clear  that  the  elements  of  P _ (0) 

nn  nn 

can  be  computed  directly  from  G(s):  that  is, 

H(0)  =  lim  s  G(8)  , 

8-»  00 

HtO)  =  lim  s[sG(s)  -  H(0)]  , 

8  -•  00 

etc. 


(3)  Compute  the  rank  n  of  F _ (0)  .  From  Theorem  4.  i, 

nn 

n  is  the  order  of  the  minimal  realizations  of  G(s). 


(4)  Let  T  be  the  submatrix  formed  from  the, first  m(n+l) 


rows  and  the  first  r(n+l)  columns  of  F 
An  explicit  minimal  realization  of  G(s) 


_ (0)  (i.e.,  T  =  F 

n  n 

is  obtained  from  T  as  in 


n+1,  n+1 


(0)). 


Theorem  4.  3. 


Vr'v.' 
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The  method  for  realizing  fixed  transfer  function  matrices 
is  quite  similar  to  those  recently  presented  by  Ho  and  Kalman  [  11  ] , 
and  by  Youla  and  Tissi  [  12  ]  while  the  present  work  was  in  progress. 
In  both  [11]  and  [  12  ] ,  a  sequence  of  matrices  equivalent  to 
FkktO)  is  generated  from  G(s)  but  the  methods  of  obtaining  a 
minimal  realization  from  these  matrices  differ  in  several  aspects 
from  that  presented  here. 


*tv 
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CHAPTER  6 


CONTROLLABILITY  AND  TIME- VARIABLE  UNILATERAL  NETWORKS* 
6.  1  Introduction 

It  has  been  seen  in  the  previous  chapters  that  the  concepts  of 
controllability  and  observability  are  valuable  characterizing  properties 
in  the  general  theory  of  linear  systems.  Although  this  viewpoint  has  not 
received  much  explicit  attention  by  circuit  theorists,  it  appears  that 
significant  problems  in  network  theory  might  be  formulated  clearly  and 
solved  effectively  through  arguments  utilizing  controllability  [  1, 2,  3,  4], 

The  field  of  time-variable  network  theory  could  especially 
benefit  from  the  introduction  of  new  techniques  since  relatively  little 
progress  has  been  made  in  this  field  by  comparison  with  the  accomplish¬ 
ments  of  fixed  network  theory.  These  considerations  motivate  a  closer 
examination  of  controllability  criteria  as  applied  to  variable-parameter 
systems  and  an  attempt  to  exploit  these  criteria  in  a  specific 
application  to  time -variable  network  theory. 

To  demonstrate  the  applicability  of  the  controllability  concept 
to  network  theory,  it  is  used  to  develop  a  class  of  networks,  composed 
of  only  linear  two-terminal  elements,  which  exhibit  unilateral  trans¬ 
mission  between  ports.  1  he  fact  that  such  networks  exist  does  not 
appear  to  be  well  known.  Although  it  is  felt  that  workers  in  the  field 

Parts  of  this  chapter  have  appeared  in  a  paper  by  the  author  and 
H.  E.  Meadows  [  8  ]  . 
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of  sampled-data  systems  may  be  aware  of  the  possibility,  to  the 
author's  knowledge  only  one  example  of  this  type  of  network  has  been 
presented  in  the  literature.  R.  J.  Mohr  [  5  ]  recently  described  a 
circuit  containing  delay  lines  and  periodic  switches  which  behaves 
essentially  as  a  three-port  circulator.  Such  unilateral  networks  have 
not  been  studied  in  any  generality,  however. 

In  this  paper  the  normal  form  of  the  network  equations  is 
used  to  show  that  unilateral  behavior  and  controllability  are  related  and 
that  the  conditions  for  controllability  developed  for  general  time-variable 
systems  can  be  used  to  characterize  classes  of  unilateral  networks. 

These  networks  must  contain  time- variable  components.  An  example  of 
a  particular  class  of  unilateral  networks  is  >resented  in  the  sub¬ 
sequent  discussion. 

6.  2  Output  Controllability 

It  will  be  shown  subsequently  that  unilateral  behavior  in  a 

network  is  intimately  related  to  the  concept  of  output  controllability 

introduced  by  Kreindler  and  Sarachik  [6].  Several  types  of  output 

controllability  may  be  defined  (as  was  done  for  slate  controllability  in 

Chapter  2)^  but  the  following  will  suffice  for  ou;  purpose  here. 

Definition  6.  1:  System  (2.  1)  is  completely  output  controllable  at  t^  , 

if  for  any  initial  state  x-,  at  t  ,  and  any  desired  final  output  y  ,  , 

U  U  d 

there  exists  a  finite  time  tj  >  t^  and  an  input  u(t)  defined  on  the 

interval  [  t,,,t,]  such  that  y(t,)  =  y,  . 

U  1  id 

The  necessary  and  sufficient  condition  for  this  type  of 
controllability  established  by  Kreindler  and  Sarachik  [  6  ]  is  stated  below. 


'■r 
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Thtorem  6.  1:  System  (2.  1)  is  completely  output  controllable  at  t^ 
if  there  exists  a  finite  time  t^  >  t^  such  that  the  rows  of  H(tj,  r  ) 
are  linearly  independent  functions  of  r  on  [  t^,  t^]  . 

In  order  that  the  concept  of  output  controllability  be 
applicable  to  th'*  study  of  time -variable  networks,  it  is  necessary  to 
obtain  a  criterion  that  does  not  depend  on  the  impulse  response  matrix. 
This  is  accomplished  below  with  the  aid  of  the  techniques  employed  in 
Chapter  3  for  state  controllability. 

Theorem  6.  2:  System  (2.  1)  is  completely  output  controllable  at  t^ 

if  for  some  t  >  t-  the  matrix  C(t)Q  (t)  has  rank  m  . 

U  c 

To  prove  this  theorem,  let  W(t,  T  )  be  the  Wronskian  matrix 
of  the  rows  of  H(t,  T  ),  that  is 

wit,T)  =  (A  (t,T) ;  A  (t,T) :  .  .  .  :a  (t,T)  ]  , 

u  .1  .  .  n- 1 

where 

^  '^,(t.T)  =  H(t,T)  . 

From  (3,  8)  it  is  clear  that 

W(t,T)  =  C(t)Mt.r)Q^(r)  . 


rygi<yyy»^ 
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Suppose  that  the  system  is  not  completely  output  controllable. 
Then,  Theorem  6.  1  implies  that  for  all  t  >  t^  ,  the  rows  of  H(t,T)  are 
linearly  dependent  functions  of  T  on  [  t^,  t  ]  .  It  follows  then  from 
Theorem  3.  1  that  W(t,T)  has  rank  <m  for  all  t  ,  and  for  all 
^q*  particular, 

W(t,t)  =  C(t)Q  (t) 
c 

has  rank  <m  for  all  t  ,  which  establishes  the  theorem. 

6.  3  Unilateral  Time-Variable  Networks 

In  the  preceding  discussion  of  controllability,  the  particular 
type  of  physical  system  represented  by  the  system  of  differential 
equations  was  not  considered.  The  present  section  is  concerned  with 
the  system  representation  of  linear  electrical  networks:  in  particular, 
the  network  property  of  unilateral  transmission  and  the  system 
property  of  controllability  are  related.  The  sufficient  conditions  for 
controllability  are  used  to  demonstrate  the  existence  of  a  class  of 
unilateral  time-variable  networks  composed  of  lumped,  linear,  two- 
tenr.inal  elements.  The  controllability  argument  ends  insight  and 
some  generality  to  this  demonstration  while  also  obviating  the 
integration  of  time- variable  differential  equations. 

The  very  well-known  concept  of  unilateral  transmission  is  a 
simple  one  v/hich  may  be  stated  in  terms  of  the  impulse  response 
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associated  with  a  pair  of  ports.  Consider  a  network  N  containing  only 
linear  two-terminal  elements,  and  let  i  and  j  be  any  two  ports  of  N. 


Denote  by  h_(t,  T}  the  voltage  [  current] ,  response  at  port  i  and 


time  t  to  a  unit  impulse  of  current  [  voltage]  at  port  j  and  time  T  . 


Then  the  network  is  a  reciprocal  network  if  h..(t,  T)  =  h, ,(t,  T)  for 

J  ^  ^3 


all  t,  r  and  all  i,  j  .  It  is  well  known  that  unless  all  of  its  elements 
are  fixed,  N  generally  may  be  nonreciprocal.  The  existence  of 
nonreciprocal  effects  suggests  the  possibility  of  completely  unilateral 


behavior  in  time-variable  networks.  If  h_(t,  t)  =  0  for  all  t  and  T  , 


and  if  for  some  t  and  T,  h..  (t,  T)  /  0  ,  N  is  unilateral  from  port  i 

3  ^ 


to  port  j  .  Networks  of  linear  two-terminal  elements  exhibiting 
unilateral  transmission  between  ports  appear  not  to  have  been 
examined  in  any  generality. 


The  main  problem  in  studying  this  type  of  unilateral 
behavior  is  that  the  networks  exhibiting  it  must  contain  time-variable 
elements.  Since  the  exact  behavior  of  a  time-variable  network  is 


difficult  to  determine,  it  seems  useful  to  develop  techniques  which 
do  not  require  knowledge  of  network  solutions.  If  the  normal  form 
(A-matrix  description)  is  used  to  write  the  network  equations,  the 
unilateral  property  may  be  phrased  in  terms  of  controllability  of 
a  linear  system  and  the  test  given  in  the  preceding  section  may  then 
be  employed  to  bypass  the  system  solution. 

The  equilibrium  equations  of  a  linear  lumped  RLC 
network  may  be  written  in  the  it  -rmal  form 


I 


•A, 


■:  -f: 
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x(t)  =  A(t)x(t)  +  B(t)u(t)  (6.1) 

as  in  (2.  1).  For  convenience,  and  with  some  loss  of  generality,  it 

will  be  assumed  that  the  input  and  output  ports  are  chosen  so  that  the 

n  vector  x(t)  represents  both  the  state  variables  (e.  g.  ,  capacitor 

voltages  for  an  RC  network)  and  the  port  outputs,  and  the  n  vector 

u(t)  represents  the  inputs  to  the  same  ports  (e.  g.  ,  current  sources 

in  shunt  with  ports  of  an  RC  network).  The  work  o£  Bryant  [  7  ]  may 

be  used  to  determine  the  matrices  A(t)  and  B(t)  explicitly. 

To  deal  with  transmission  between  pairs  of  ports  of  the 

network,  it  is  convenient  to  recast  (6.  1)  in  a  slightly  modified  form. 

th 

Let  e^  represent  the  unit  column  vector  with  a  one  in  the  i  position 
and  zeros  elsewhere,  and  e!  the  corresponding  unit  row  vector;  then 
in  order  to  consider  inputs  applied  only  at  port  i  and  outputs  only  at 
port  j  ,  (6.  1)  may  be  replaced  by  the  system  F..  described  by 


x(t)  =  A(t)x(t)  +  B(t)e^u(t) 

y(t)  =  x(t)  (6.2) 


Theorem  6.  3:  The  network  with  normal  form  equations  (6.  2)  is 
unilateral  from  i  to  j  if  and  only  if  system  F_  is  nowhere 
controllable,  and  system  is  completely  output  controllable  at  some 
time  tjj. 
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By  "nowhere  controllable"  we  mean  that  the  system  is  not 
completely  output  controllable  at  any  time.  The  theorem  is  clear 
since  h  (t.^ )  is  simply  the  impulse  response  of  F..  ;  that  is, 

h.^(t.T)  =  e:f(t,T)B(r)e.  . 

To  establish  that  networks  satisfying  the  criteria  of  Theorem 

6.  3  do  exist,  a  class  of  RC  unilateral  networks  will  be  derived. 

- 1 

Provided  that  C  (t)  and  C(t)  exist  for  all  t  ,  the  equilibrium 
equations  for  an  RC  network  can  be  written  as 

V(t)  +  C“Nt)[6(t)  +  G(t)]  V(t)  =  C■^t)J(t),  (6.3) 

where  V(t)  represents  a  set  of  voltage  variables  pertinent  to  a  tree  in 
a  capacitive  subgraph  of  the  network,  and  C(t)  and  G  (t)  are  the  n- 
port  capacitance  and  conductance  matrices  corresponding  to  the  voltage 
variables.  The  current  vector  J(t)  represents  a  set  of  current 
sources  in  parallel  with  the  ports  determined  by  the  voltage  variables. 
Thus,  the  order  of  the  system  (6.  3),  is  the  same  as  the  order  of  com¬ 
plexity  of  the  network,  and  the  elements  of  V(t)  and  J(t)  represent 
output  voltage  and  input  currents,  respectively.  It  will  now  be  shown 
that  a  network  with  system  equation  (6.  3)  may  be  unilateral.  Define 
S(t)  =  C  ^(t);  then  system  F. .  corresponding  to  (6.3)  is 
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V(t)  =  -S(t)[  C(t)  +  G(t)]  V(t)  +  S(t)c.j.(t)  (6.  4) 

J  J 


v.(t)  =  e!V(t)  . 


A  necessary  condition  for  F. .  to  be  nowhere  controllable  is  that  the 

ij 

ij  element  of  S(t)  must  vanish  identically.  If  the  contrary  were  true 
then  h.j(t,r)  =  0  would  imply  that,  for  all  t  and  for  all  r  such  that 
s..(T)  i  0  . 


1 

s..(r) 


£ 

k=l 


(6.5) 


But, 


(T,r) 


0  ,  k  i 

1  ,  k  =  i  ; 


therefore,  (6.5)  is  contradictory  at  t  =  T  and  must  be  zero. 

^  _  Consider  now  the  case  n  =  3  ,  and  let  i  =  3  and  j  =  1  . 

Since  the  port  matrices  C(t)  and  G(t)  are  symmetric. 


i 

Si3(t)  =  8^j(t)  =  0  . 


With  this  constraint,  the  condition  for  Fjj  to  be  output  controllable  at 
some  tp  (Theorem  6.  2)  is  that  the  matrix 
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ej'[  0  :  -A(t)S(t)ej  !  {  A^(t)S(t)  -  2A(t)S(t)  -  A(t)S{t)}e3]  (6.  6) 

have  rank  1  for  lome  t,  >  .  k4oreover,  F, ,  is  nowhere 

10  31 

controllability  if  and  only  if 


hjjU.T)  =  ‘?>3i(t.T)8jj(T)  + 

To  simplify  the  presentation  without  great  loss  of  generality,  let 

812(0  =  82j(t)  =  0  : 


then  (6.7)  becomes 


''3l(t,T)8ji(r)  =  0. 


(6.8) 


The  A(t)  matrix  for  system  (6.  4)  can  now  be  written  as 
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*13 


*^23 


^  0  . 


(6.  10c) 


*^22^^  23  ^ 


22*^23  • 


(6.  lOd) 


An  illustrative  network  is  provided  by  the  following 

example. 

Example  6,  1;  Consider  the  circuit  of  Fig.  6.  1,  whose  port  matrices 
are 


G  = 


1 

4 

4  \  cost 

5  +  cost 


4  +  cost 
sint 
0 


/  ^  cost , 
5  +  cost 

0 

1 


C  = 


0 

5  +  cost 
-1 


Examination  of  G  and  C  shows  that  the  constraints  in  (6.  10)  are 
satisfied,  so  that  the  network  transmits  from  port  3  to  port  1,  but  not 
in  the  opposite  direction.  It  may  be  observed  that  the  constraints  in 
(6.  10)  are  independent  of  Cjj,  g^j,  c^j,  and  g^^.  Therefor^  ports  1 
and  3  could  be  padded  arbitrarily  while  preserving  the  unilateral  trans- 
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mission  property,  as  would  be  expected  intuitively.  Several  other 
examples  of  RC  unilateral  networks  are  given  in  f  8  ]  . 

The  network  of  Example  6.  1  contains  one  negative  element, 
but  as  will  be  shown  below,  the  circuit  behaves  in  a  stable  manner  under 
all  excitations.  This  analysis  will  also  contribute  further  physical 
insight  into  the  unilateral  nature  of  the  network. 

Consider  the  network  of  Example  6.  1  with  the  elements  in 
parallel  with  ports  1  and  3  removed  as  shown  in  Fig.  6.  2.  This 
simplified  circuit  may  also  be  represented  as  in  Fig.  6.  3. 

It  may  be  readily  verified  that  if  the  network  is  excited  by 

jj  ,  then 


^2  =  + 


(6.  11) 


and 


V 


4 


-  b(t)jj  , 


(6.  12) 


where 


b(t) 


1 

4  cost 


It  is  clear  that  for  any  input  ,  the  output  voltage  at  port  3  is 


A 


.  ,  ... 


iU 


identically  zero  for  zero  initial  conditions,  since 


Vj  =  =  0. 


(6.  13) 


If  the  network  is  excited  by  however,  this  type  of  voltage 
cancellation  does  not  take  place  at  port  1  since  is  a  nonreciprocal 
network,  as  will  now  be  shown.  The  response  at  port  1  of  is 
given  by 


'4  =  •'’4*j3 


Vj  .  -bv^ 


(6.  14) 


v,  =  -('  -b ’"S  ■  '=^3  • 


(6.  15) 


while  that  at  port  2  of  N,  is  the  same  as  (6.  11)  with  j  replaced  by 

b  1 

.  The  nonreciprocal  nature  of  is  evidenced  by  the  term  b/b  . 
If  b  were  constant,  (6.  15)  would  be  identical  with  (6.  12).  Thus, 


v  .  +  v_  /  0  . 

4  5 
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The  impulse  response  in  the  direction  of  non- zero  transmission  is 
easily  found,  and  is  given  by 

hj2(t,r)  =  e"^^‘’’^b(r)  -  b(t)e‘^^‘’‘^  .  (6.16) 

It  is  clear  from  (6.  16)  that  the  above  network  is  stable  in  the  sense 
that  every  bounded  input  j^lt)  will  produce  a  bounded  output  Vj(t) 

[  9  ]  ,  since  for  all  t 

7  <  b(t)  <  1  . 

Equations  (6.  11)  and  (6.  12)  are  similarly  stable. 

The  fact  that  unilateral  transmission  is  possible  in  networks 
containing  only  time- variable  RC  elements  indicates  that  many  practical 
devices  might  be  devised  employing  the  nonreciprocal  effects.  Current 
investigations  not  yet  completed  have  yielded  configurations  whose 
terminal  behavior  is  equivalent  to  that  of  such  bas  c  circuits  as  gyrators, 
nullators,  and  norators  [  13  ]  .  Although  these  realizations  can  be  shown 
to  exist  mathematically,  many  practical  problems  remain  before  they 
can  be  implemented  physically.  Chief  among  these  are  the  feasibility 
of  the  required  element  variations,  and  the  sensitivity  of  the  unilateral 
transmission  to  perturbations  in  the  elements. 

Another  possible  area  of  application  might  be  ir  the  design 
of  parametric  amplifiers.  In  recent  years  several  papers  have  dealt 
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with  unilateral  parametric  amplifiers  not  requiring  unilateral  devices 
[lO,  11,  12],  These  amplifiers  behave  unilaterally  only  for  fixed 
excitation  frequency,  however,  and  when  viewed  as  linear  circuits  with 

i  their  sources  removed  they  do  not  have  the  unilateral  property  as 

defined  here. 


i 
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Chapter  7 
CONCLUSION 

In  this  thesis,  a  new  characterization  of  several 
fundamental  linear  system  concepts  has  been  presented  which 
employs  only  explicitly  known  system  parameters  and  does  not 
require  the  solution  of  time-variable  differential  equations.  With 
this  characterization,  it  was  possible  to  provide  a  unified  approach  to 
the  representation  and  synthesis  of  a  broad  class  of  variable 
parameter  systems,  and  to  derive  an  interesting  class  of  time- 
variable  networks.  A  detailed  summary  of  the  major  results  obtained 
is  given  below. 

7.  1  Summary  of  Results 

It  has  been  shown  that  the  problem  of  determining  the 
degree  of  controllability  or  observability  of  a  time-variable  linear 
system  in  terms  of  its  coefficient  matrix  description  is  essentially 
equivalent  to  that  of  determining  the  degree  of  independence  of  a  set 
of  functions  by  a  Wronskian  type  test.  By  developing  such  a  test 
for  vector  functions,  precise  criteria  for  controllability  and 
observability  have  been  obtained  which  do  not  require  the  solution  of 
time- variable  differential  equations.  These  criteria  are  based  on 
the  rank  of  the  controllability  and  observability  matrices,  newly 
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defined  quantities  formed  from  the  system  coefficient  matrices  and 
a  finite  number  of  their  derivatives.  The  controllability  and 
observability  matrices  were  also  shown  to  be  useful  in  a  variety  of 
other  system  analysis  problems.  In  particular,  a  transformational 
property  of  these  matrices  was  employed  which  enabled  the 
formulation  of  general  criteria  for  system  reducibility  and  led  to 
powerful  methods  for  reducing  noncontrollable  or  nonobservable 
systems  to  lower  order.  Also  obtained  were  new  methods  for 
characterizing  and  generating  equivalent  system  representations. 
Criteria  for  equivalence,  zero>state  equivalence,  and  zero- state 
time- invariance  of  time-variable  linear  systems  were  among  the 
results  of  this  development.  Several  useful  canonical  forms  for  time- 
variable  linear  systems  were  also  investigated  in  this  context  and 
criteria  for  their  existence  and  method  for  their  construction  were 
given. 

Based  on  the  theory  of  equivalent  systems  and  the  criteria 
for  controllability  and  observability,  a  new  approach  to  the  synthesis  of 
nonstationary  impulse  response  matrices  was  developed.  This  method 
of  synthesis,  which  does  not  require  an  a  priori  assumption  of 
separability,  provides  a  systematic  procedure  for  realizing  a  wide 
class  of  system  impulse  responses.  The  specialization  of  this  pro¬ 
cedure  to  fixed  systems  is  of  independent  interest,  since  it  has 
several  advantages  in  comparison  with  existing  techniques. 
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By  relating  the  concepts  of  controllability  and  unilateral 
transmission  it  was  possible  to  use  the  controllability  criteria  to 
demonstrate  the  existence  of  a  class  of  unilateral  networks  composed 
solely  of  two-terminal  RC  (time- variable)  components.  A  stable 
example  of  such  a  network  was  presented  and  possible  applications 
discussed. 

7.  2  Problems  for  Future  Research 

(1)  It  is  well  known  [1,2]  that  the  optimvun  filter  for 
determining  the  state  variables  of  a  system  from  observations  of  its 
output  requires  differentiators  when  some  or  all  of  the  output 
measurements  are  assumed  to  be  noise  free.  The  precise  conditions 
for  existence  of  such  a  filter  have  not  been  given,  however,  nor  the 
details  of  its  construction.  In  the  total  absence  of  noise,  equation 
(4.  8)  essentially  provides  the  form  of  the  filter,  and  it  is  clear  that 
uniform  observability  is  the  required  condition  for  existence.  A 
detailed  examination  of  the  observability  matrix  should  provide  a 
method  for  constructing  such  filters  and  determining  whether  they 
exist  \inder  the  more  general  noise  constraints  considered  by  Bryson 
and  Johansen  [  2  ]  . 

(2)  The  solution  of  many  fixed  control  system  problems  is 
greatly  simplified  by  first  transforming  a  given  system  to  an 
appropriate  canonical  form  f  3-6  j  .  In  fact,  the  conditions  for 
existence  of  solutions  to  most  of  these  problems  coincide  with  the 
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conditions  for  existence  of  the  canonical  equivalents.  It  is  expected, 
therefore,  that  the  theory  and  methods  of  constructing  canonical 
forms  developed  in  Chapter  4  will  prove  to  be  valuable  in  extending 
previously  restricted  results  to  time-variable  systems.  For 
example,  it  can  be  shown  with  the  aid  of  form  (4.  33)  that  if  a  fixed 
single-input  system  is  completely  controllable,  any  desired  pole 
pattern  may  be  achieved  by  feeding  back  a  linear  combination  of  the 
state  variables  of  the  system.  While  the  concept  of  a  "pole"  is 
somewhat  nebulus  for  time-variable  systems,  it  is  clear  that  if  all  the 
coefficients  a^  in  (4.  33)  can  be  adjusted  independently  by  feedback, 
the  system  behavior  can  be  modified  as  freely  as  in  the  fixed  case. 

If  a  system  is  uniformly  controllable,  transformation  to  form  (4.  33) 
is  possible  and  the  coefficients  of  the  required  feedback  vector  can 
be  found  precisely  as  in  f  3-4  ]  . 

(3)  It  would  be  desirable  to  extend  the  synthesis  procedure 

N.. 

of  Chapter  5  to  more  general  impulse  response  matrices.  While  it 
seems  unlikely  that  all  linear  systems  can  be  realized  by  this  method 
it  should  be  possible  to  include  realization  of  impulse  response  matrices 
of  total  controllable  and  observ'able  systems.  In  addition,  the  use  of 
the  present  method  as  part  of  an  approximation  procedure  for 
realizing  numerically  tabulated  response  data  (c.8  would  usually  be 
specified  in  practice)  should  be  investigated.  A  detailed  evaluation 
of  the  various  method  for  synthesizing  fixed  responses  should  also  be 
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undertaken  to  determine  the  one  most  suitable  for  machine 
computation. 

(4)  As  indicated  in  Chapter  6,  there  are  many  p.  oblems 

remaining  before  practical  unilateral,  and  other  nonreciprocal 
networks  can  be  realized  with  practical  time-variable  components. 
If  these  problems  can  be  overcome,  however,  a  new  method  for 


building  the  non  reciprocal  elements  widely  used  in  network  theory 
would  Le  available. 
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APPENDIX 


PROPERTIES  OF  ANALYTIC  FUNCTIONS 

The  following  well  known  properties  of  analytic  functions  are 
of  importance  here. 

Property  A.  1 ;  If  two  analytic  functions  are  equal  et  more  than  a  finite 
number  of  points  on  a  finite  interval,  they  are  equal  everywhere  on  the 
interval.  In  particular,  an  analytic  function  that  is  zero  at  more  than 
a  finite  number  of  points  on  a  finite  interval  is  identically  zero  on  the 
interval. 

1 

Property  A.  2;  The  sum  and  product  of  two  analytic  functions  is  an 
analytic  function. 

Property  A.  3;  If  an  analytic  function  is  non  zero  on  an  interval,  its 
inverse  is  analytic  on  the  interval. 

Property  A.  4;  The  derivative  of  an  analytic  function  is  an  analytic 
function. 

The  above  properties  of  scalar  analytic  functions  will  nov 
be  utilized  to  derive  several  properties  of  matrices  of  analytic  functions 
(analytic  matrices). 

Property  A.  5:  The  determinant  of  any  square  submatrix  of  an  analytic 
matrix  is  an  analytic  function. 

Proof:  Follows  directly  from  Property  A.  2. 
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Property  A.  6;  The  rank  of  an  analytic  matrix  is  constant  save  possibly 
at  a  finite  number  of  points  on  any  finite  interval. 

Proof;  Let  f(t)  be  the  determinant  of  the  largest  square  submatrix 
nonsingular  for  some  point  on  a  given  finite  interval.  It  follows  from 
Properties  A.  1  and  A.  5  that  f(t)  must  be  nonzero  except  possibly 
at  a  finite  number  of  points  on  the  interval. 

Property  A.  7;  If  an  analytic  matrix  has  rank  q  save  possibly  at  a 
finite  tiximber  of  points  on  a  finite  interval  then  a  square  submatrix  of 
order  q  also  has  this  property. 

Proof;  Follows  directly  from  Property  A.  6. 

It  is  now  possible  to  prove  Corollary  3.  2  of  Chapter  3. 

First  note  that  Property  A.  4  implies  that  W(t)  is  an  analytic  matrix 

if  the  row  vector  functions  0, ,  6  , ....  6  a^re  analytic. 

1  c  n 

If  W(t)  has  rank  n  for  some  t  €  [  t^,  tj]  then  Theorem  3.1  implies 
that  the  functions  are  independent  on  the  interval.  Suppose  then  that 
W(t)  has  rank  leas  than  n  for  all  t  e  [  which  case  it  must 

have  an  essentially  constant  rank  q  <  n  on  the  interval  (Property  A.  6). 
Let  U  be  a  SLlinterval  of  [  t^,  tj]  over  which  q  rows  of  W(t) 
have  rank  q  everywhere  (that  such  a  subinterval  exists,  follows 
from  Property  A.  7),  From  Theorem  3.2,  there  exists  a  constant 
nonsingular  matrix  such  that  on  U 


rw(t) 

TWt^)  --  i  ,  (A.  1) 

L  0  J 


where  W(t)  has  q  rows  and  rank  q  .  Since  TW(t)  is  analytic 
on  the  entire  interval  [  t^]  the  relationship  (A.  1)  must  also  hold 

on  [iQjij]  (Property  A.  1).  Thus,  as  in  Theorem  3.2  the  6^  are 
dependent  over  [  ^q,  tj]  while  exactly  q  of  the  0^  are  linearly 
independent  on  the  interval  (i.  e. ,  those  corresponding  to  the  rows 
of  W(t)  having  lank  q  on  U  ). 
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